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ABSTRACT 

We examine the properties of the diffuse hard X-ray emission in the classic starburst galaxy M82. We 
use new Chandra ACIS-S observations in combination with reprocessed archival Chandra ACIS-I and XMM- 
Newton observations. We find E ^6.1 keV Fe Hea emission (from highly ionized iron) is present in the central 
\r\ < 200 pc, |z| < 100 pc of M82 in all datasets at high statistical significance, in addition to a possibly non- 
thermal X-ray continuum and marginally significant E = 6A keV Fe Ka line emission (from weakly ionized 
iron). No statistically significant Fe emission is found in the summed X-ray spectra of the point-like X-ray 
sources or the ULX in the two epochs of Chandra observation. The total nuclear region iron line fluxes in the 
2004 April 21 XMM-Newton observation are consistent with those of the Chandra-derived diffuse component, 
but in the 2001 May 6 XMM-Newton observation they are significantly higher and also both E=6.4 and E=6.9 
keV iron lines are detected. We attribute the excess iron line emission to the Ultra-Luminous X-ray source 
in its high state. In general the iron K-shell luminosity of M82 is dominated by the diffuse component. The 
total X-ray luminosity of the diffuse hard X-ray emission (corrected for emission by unresolved low luminosity 
compact objects) is Lx,2-8kev ^ 4.4 x 10^^ erg s~' in the £ = 2 - 8 keV energy band, and the 6.7 keV iron line 
luminosity is Lx.e.vkev ~(1.1 - 1-7) x lO-** erg s"' . The diffuse hard X-ray continuum is best fit with power law 
models of photon index T ^ 2.5 - 3.0, although thermal bremsstrahlung model with kT ^ 3 - 4 keV are also 
acceptable. A simple interpretation of the hard diffuse continuum is that it is the bremsstrahlung associated 
with the 6.7 keV iron line emission. This interpretation is problematic as it would imply an unrealistically 
low gas-phase iron abundance Zpe ~ QAZp^ Q, much lower than the super-Solar abundance expected for su- 
pernova ejecta. We explore the possibility of non-equilibrium ionization, but find that this can not reconcile 
the continuum and iron line emission. Nor can Inverse Compton X-ray emission be the dominant source of 
the diffuse continuum as theory predicts flatter power law slopes and an order of magnitude less continuum 
flux than is observed. The 6.7 keV iron line luminosity is consistent with that expected from the previously 
unobserved metal-enriched merged supernova ejecta that is thought to drive the larger-scale galactic superwind. 
The iron line luminosity implies a thermal pressure within the starburst region of P/k ^2x lO^Kcm"^ which 
is consistent with independent observational estimates of the starburst region pressure. 

Subject headings: ISM: jets and outflows — ISM: bubbles — galaxies: individual : NGC 3034 (M82) — 
galaxies: halos — galaxies: starburst — X-rays: galaxies 



1. EVITRODUCTION 

Starburst galaxies - objects with intense recent star forma- 
tion - have galaxy-sized outflows with velocities in excess of 
several hundred to a thousand kilometers per second. These 
superwinds are common in the local Universe, occurring in 
nea rly all galaxies classified as undergoing sta rburst activ- 
ity (iLehnert & Heckmanlll995l: lHeckmanlll998h . and appeai- 
ubiquitous among the Lyman break galaxies at redshift ~ 3 
where they are blowing ^100 kpc sized hole s in the IGM 
(e.g. Pettini et alJl200iriAdelberger et al.ll2003h . 

Superwinds are believed to be driven by the thermal and 
ram pressure of an initially very hot (T ^ 10^ K), high pres- 
sure (P/k ^ lO^Kcm"^ ) and low density wind, itself created 
from the merged remnants of very large numbers of core- 
collapse supernovae (SNe), and to a lesser extent the stel- 
lar winds from the massive stars, that oc cur over the ~ 100 
Myr d uration of a typical starburst event (IChevalier & Cleg^ 
Il985h . Direct observational detection of gas this hot through 
X-ray spectroscopy or imaging has only recently become pos- 
sible because of the high spatial resolution of Chandra X- 
ray Observatory, as starburst regions are also home to large 
numbers of X-ray-luminous accreting neutron stars and black 
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holes that must be resolved out in order to see any diffuse hard 
X-ray emission. 

All known superwinds were detected through observa- 
tions of ambient disk or halo gas that has been swept- 
up by the hotter but more tenuous wind of merged SN 
ejecta^, which has itself escaped direct observational de- 
tection until recently. The majority of nearby super- 
winds have been discovered using optical imaging and spec- 
troscopy to identify outflow in the warm ionized medium 
(WIM, r '--^ 10"* K), in many cases directly imaging bipo- 
lar structures aligned with the host galaxy minor axis 
that have kinematics indicative of outflow at velocities of 
v^^iM = a few xlOO to lOOOkms"' (e.g. lAxon & Tavloii 
1978; Heckman, Armus, & Miley 1987; Bland & Tully l988|; 
Heckman, Arm us. & Mile v 1990; Lehnert & Heckman 19951 
1996). McCar thv et all (fl 987), Heckman etal. (1990), and 
ILehnert & HeckmanI (11996 ) use this entrained low-filling 
factor gas as a tracer of the pressure in the inner kilo- 
parsec or so of many nearest brightest superwinds, and 
find the very high central pressures and radial variation in 

^ It is useful draw a conceptual distinction between the material that actu- 
ally drive superwinds, and the many other phases of swept-up and entrained 
gas that are part of superwinds. We shall define the wind-fluid as the merged 
SN ejecta (possibly mass-loaded by ambient gas) whose high thermal and/or 
ram pressure drives the superwind. 
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pressure expected from the IChevalier & Clegj d 1985b wind 
model for the thermalized SN ejecta within the vicinity 
of the starburst region. Superwinds from galaxies at high 
redshift are recognized by virtue of blue-shifted interstel- 
lar absorption lines from warm-neutral and wa rm-ionized 
gas species (Pettini et alJ l2000t iFrye et all 12002) . absorp- 
tion features very similar to those seen in local starburst 
galaxies with superwind s ('Phinips"1993'; 'Kunthet al.l ll998l: 
iGonzalez Delgado et al.] [T998; Heckman et al. 2000). Nearby 
superwinds have also been extensive l y studied in soft X-ray 
emission-' (e.g. 

iDahlem et al.ll l998t Istrickland et al.||2004|) 
which probes emission from hot gas with temperatures in 
the range T ~ 10'' to lO^K, fol lowing their initial detection 
by Einstein X-ray observatory jWatson. Stanger. & Griffithsl 
[1984; Fabbiano & Trinchieril 11984 . Indeed, observations 
have demonstrated that all phases of the ISM found in normal 
late type galaxies are also in corporated into starburst-driven 
superwinds (see the review of Dahlem|[T997l) . 

However, prior to the launch of the Chandra X-ray Obser- 
vatory (CXO) in 1999, there had been no believable detec- 
tion of thermal X-ray emission from the T ^ a. few xlO^ 
to 10** K merged SN ejecta whose thermal pressure initially 
drives the superwind in any starburst galaxy. This very 
hot plasma should be a source of thermal bremsstrahlung 
in the hard X-ray energy band, as well as strong £ « 6.7 
and (possibly) 6.9 keV emission from highly-ionized Helium- 
like and Hydrogen-like iron. Hard X-ray emission had 
been previously been detected from nearby starburst galax- 
ies such as M82 and NGC 253 w ith X-ray observatories such 
as EXOSAT (Schaafet al."1989'), G/n^a (Ts uru et all 1 19901) . 
ASCA (Ptak et al. 1997; Tsuru et al. 1997) and BeppoSAX 
dPersic et al.lll998l;ICappi et al.lll999l) . However, the poor- to- 
nonexistent spatial resolution of these observations prevented 
a clean assessment of the relative contribution to the hard X- 
ray emission from compact accreting objects and from any 
truly diffuse component (the half power enclosed angular di- 
ameter of the point spread function (PSF) of BeppoSAX and 
ASCA at £■ ^ 6 keV were rso ^ 2' and ^ 3' respectively. This 
was the highest spatial resolution available for any hard X-ray 
observatory before the launch of Chandra. Assuming a dis- 
tance of 2.6 Mpc to NGC 253 and 3.6 Mpc to M82, an angular 
size of 2' corresponds to a physical size of 1 .5 kpc for NGC 
253 and2.1kpcforM82). 

Nevertheless it was clear that the majority of the hard X- 
ray emission from these galaxies came from compact objects. 
The observed total hard X-ray luminosities of these galaxies 
were > 1 .5 orders of magnitude higher t han expected from the 
the ve ry hot gas in a superwind model ( Strickland & StevensI 
120001) . Both the ASCA and BeppoSAX observations of M82 
and NGC 253 also discovered emission possibl y due to high 
ionized iron (Ptak et al. 1997; Tsuru et al. 19971 iPersic et all 
[1998'; Cap pi et al.lll999D . although this emission was of low 
equivalent width. The weakness of the Fe K-shell emission 
implied sub-Solar Fe abundances {Zfe < 0.3Z/r(, 0), a pecu- 
liar result given that super-Solar abundances {Zf^ ^ SZpe.Q) 
are expected from pure core-collapse SN ejecta. Alternatively 
the hard X-ray emission could be dominated by non-thermal 
continuum emission, possibly from X-ray binaries, thus re- 
ducing the equivalent width of the Fe line from the diffuse hot 

^ We will use the following definitions in this paper: Soft X-ray denotes 
photon energies in the range £ ~ 0.1 - 2.0 keV. Hard X-ray denotes photon 
energies E > 2keV, and within the range 2 < E < lOkeV unless it is stated 
otherwise. 



gas. Furthermore, some Galactic X-ray binaries do sh ow Fe 
K-shell emission (e.g. see IVrtilek. Soker. & Ravmondiri 993). 
so the Fe emission observed in M82 and NGC 253 might 
come from compact objects rather than any diffuse gas. De- 
spite the low spatial resolution of ASCA, the hard X-ray emis- 
sion in NGC 2 53 was dist r ibuted so broadly over the disk of 
the galaxy that iPtak et all ( Il997h were able to resolve mul- 
tiple individual hard X-ray sources. These were also di- 
rectly associated with point-like X-ray sources seen in soft 
X-ray observations with the ROSAT High Resolution Imager 
(HRI). The stron g temporal va riability i n the hard X-ray lu- 
mino sity of M82 (iPtak & Griffit hs 1999; Matsumoto & Tsurul 
1 19991) demonstrated that one or more accreting compact ob- 
jects dominated the hard X-ray emission at some, and possibly 
all, epochs. 

The on-axis half-enclosed-power PSF diameter of Ad- 
vanced CCD Imaging Spectrometer (ACIS) on the Chan- 
dra X-ray Observatory is 0"8, a factor of 150 smaller than 
BeppoSAX, and also considerably better than the equiva- 
lent value of 13 - 15" for the MOS and PN detectors on 
XMM-Newton. The angular diameter of the starburst region 
in M82 is - 45" (see e.g. Golla, Allen, & Kronberg 1996), 
while that of NGC 253 is < 15" (Ulvesta d & Antonuccik 19971: 
iForbes et al.ll2000l) . It is clear that Chandra is the first, and 
only. X-ray instrument with spatial resolution high enough 
to resolve the starburst regions in even the nearest starburst 
galaxies such as M82 and NGC 253. 

Using the spectral-imaging CC Ds of the ACIS-I camera 
of Chandra, ([Griffiths et al.l l2000h detected diffuse hard X- 
ray emission coincident with the nuclear starburst region of 
M82. They interpreted the spectrum of this diffuse emission 
at thermal bremsstrahlung from a hot plasma of temperature 
r ~ 4 X lO^K. This was motivated by the detection of excess 
emission in the energy range £ = 5.9 to 6.9 keV which they at- 
tributed to unresolved line emission from highly-ionized iron. 
This feature that would not be expected if the diffuse X-ray 
emission where non-thermal in origin, for example inverse 
Compton (IC) scattering of IR photons off relativistic elec- 
trons. This work was based on ACIS-I observations taken 
in late 1999, soon after the cosmic -ray-induced damage to 
front-illuminated ACIS-I CCDs had been discovered. As a 
result, the spectral calibration of these observations was espe- 
cially difficult. In particular, the radiation-induced reduction 
in spectral resolution and absolute energy scale calibration of 
the instrument made it difficult to be sure that Fe emission 
feature was indeed real. 

ACIS-S observations also uncovered diffuse hard X-ray 
emis sion in the compac t nuclear starburst region of NGC 
253 dWeaver et al.ll2002h . The emission is strongest within 
the central 5" of the galaxy, and showed a strong hard X- 
ray emission line spectrum with features attributable to Si, 
S, Ca and Fe (the back-illuminated chips on ACIS-S suf- 
fered less radiation damage than the front-illuminated chips, 
and hence have higher spectral resolution). This spectrum is 
more reminiscent of a highly photoionized plasma than a col- 
lisionally ionized plasma. As there is independent evidence 
for a low lu minosity AGN (LLAGN) in the cen ter of NGC 
253 (lUlvestad & Anton ucci[ [l997l) . [Weaver et all jlOOl) sug- 
gest that the diffuse hard X-ray emission from NGC 253 is 
dominated by gas photoionized by the LLAGN. 

Thus the nature of this diffuse hard X-ray emission in M82, 
and its possible relationship to the very hot merged SN ejecta 
thought to drive M82's superwind, has not yet been resolved 
and is the focus of this paper. We emphasize that M82 is the 
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best case for exploring the properties of diffuse very hot gas 
created by starburst activity, as the diffuse hard X-ray emis- 
sion is brighter and has a larger angular extent than in NGC 
253. Furthermore, there is no known LLAGN in M82 to com- 
plicate matters. 

In this paper we make use of the longest available Chan- 
dra ACIS and XMM-Newton observations and the most re- 
cent instrumental calibrations to further investigate the prop- 
erties of the diffuse hard X-ray emission in this archetypal 
starburst galaxy. With Chandra we spatially separate the dif- 
fuse hard X-ray emission from resolved point source emis- 
sion, and use the XMM-Newton observations provide the best 
signal-to-noise spectra of the nuclear starburst region of the 
galaxy. These observations also sample four different epochs, 
allowing us to separate the time-constant diffuse component 
of the hard X-ray iron line emission from a one-time contribu- 
tion by the Ultra-Luminous X-ray source found in the shorter 
XMM-Newton observation. 

2. OBSERVATIONS AND DATA REDUCTION 

We obtained a 18 ks-long observation of M82 on 2002 June 
18 (1061 days since the launch of Chandra), with the cen- 
tral few kpc of the galaxy on the back-illuminated S3 chips of 
the Chandra ACIS-S instrument (observation ObsID 2933). 
This ACIS-S observation has higher spectr al resolution tha n 
the ACIS-I observations presented in Griffit hs et al.l 120001). 
In addition this new observation, we have reanalyzed two 
archival ACIS-I observations of M82, both taken on 1999 
September 20 (55 days since launch, ObsIDs 1302 and 361, 
PI Murray and Garmire respectively), making use of the latest 
spectral calibrations, and corrections for the radiation-damage 
induced charge transfer inefficiency (CTI) in the AC IS-I chips 
that were unavailable at the time of iGriffiths et al.l study. The 
longer of these two ACIS-I observations (O bsID 361, ap- 
proxim ately 33 ks in length) was analyzed by IGriffiths et alj 
(l2000h . 

In addition to reanalyzing this dataset, we have also com- 
bined the ObsID 1302 and ObsID 361 data into a single longer 
dataset. Although image-based analysis of multiple merged 
Chandra datasets is common in the literature, it is normally 
inadvisable to perform spectral analysis on such merged data. 
In this particular case this is not a concern as the observa- 
tions were performed concurrently, and the pointing direc- 
tions and telescope roll of the observations are essentially 
identical. As photons within the nuclear region of interest fell 
on the same pixels of the same chips in both observations, the 
chip-location-specific instrument responses are the same and 
we can proceed with a spectral analysis on the merged data. 

For convenience we shall refer to the different datasets ana- 
lyzed here as the ACIS-S (obsID 2933), 33 ks ACIS-I (ObsID 
361) and merged ACIS-I (ObsID 361 + ObsID 1302) obser- 
vations. 

The Chandra data was reduc ed and analyzed in la rgely 
the same ways as described in [Strickland et al.l (|2004|) . al- 
though using more recent versions of the software and im- 
proved spectral and energy calibrations. A few details of the 
data reduction of specific relevance to this pa per are discussed 
here, in particular the differences from the the lStrickland et al] 
(|2004) data reduction. 

We also reduced and analyzed the data for the two longest 
XMM-Newton observations of M82, obtaining the data from 
the public archive. The EPIC instruments on XMM-Newton 
do not have high enough spatial resolution to separate the dif- 
fuse and point source related X-ray emission in the nucleus 



of M82. We make use of EPIC's higher effective collecting 
area to place high signal-to-noise constraints on the total K- 
shell iron line emission from M82 at different epochs from the 
Chandra observations. 

The Chandra and XMM-Newton observations are summa- 
rized in Table [U 

We adopt a distance of 3.6 Mpc to M82. This is based on 
the Cepheid distance of 3 .63 ± 0.34 M pc to its close neigh- 
bor galaxy M8 1 ( Freedma n et al. Il994t) . to which M82 is con- 
nected by Hi tidal tails (Cottrell 197/). This distance is con- 
sistent with detection of the tip of the red gia nt branch in M82, 
which indicates a distance of 3.9 ± 0.4 Mpc (ISakai & Madord 
Il999l) . At our adopted distance an angular size of 1" cor- 
responds to a physical size of 17.5 pc. We also adopt 
a posi tion angle of 73° fo r the major axis of the galaxy 
( Achte rmann & Lacy! 1 19951) and the location of the dynami- 
cal center of th e galaxy at a = 09''55'"5 1?9, 5 = +69°40'47."1 
(J2000.0) from lWeliachew et alJ (fl984l) . 

2.1. Chandra data 
2.1.1. CTI and gain correction 

The data was reduced using CIAO (version 3.2). The lat- 
est calibrations from version 3.0 of the Chandra calibration 
database (CALDB). Note that this version of CIAO can cor- 
rect for the radiation-induced CTI (manifest as an appar- 
ent decrease photon event energy, and spectral resolution, 
with increasing distance from the CCD readout) of the front- 
illuminated ACIS CCDs, but not the back-illuminated CCDs, 
for observations with the focal plane instruments at -120° 
C. However, the effect of CTI in the ACIS-S observations 
of M82 is mitigated by the placement of the nuclear region 
within ~ 1/3 of a chip- width from the readout edge of the S3 
chip, and the relative lack of radiation damage to the back- 
illuminated chips SI and S3. The changes in apparent pho- 
ton energy due to ongoing increases in CTI in the the back- 
illuminated CCDs is largely corrected for in CIAO 3.2 by the 
incorporation of time-dependent gain corrections'*. 

CIAO 3.2 can not CTI-correct the earlier ACIS-I observa- 
tions of M82, which were taken with a focal plane tempera- 
ture of-1 10° C. In these observations the nucleus of M82 is at 
the far edge of the 13 chip from the readout edge, maximizing 
the negative effects of CTI. We used the CTI-correction soft- 
ware (version 1.45) dev eloped by the P enn State University 
(PSU) ACIS group (see lTownsley et al.l l2000) to CTI-correct 
the M82 ACIS-I data, before reducing the data using CIAO 
as normal. This corrects for the position-dependent energy 
degradation, including the secular increase in CTI discussed 
above, but does not correct for decrease in spectral energy 
resolution also associated with CCD CTI. For this reason the 
CTI-corrected ACIS-I observations of M82 have lower spec- 
tral resolution than the ACIS-S observations. 

The standard Chandra data pipe-line introduces an artificial 
0."5 Gaussian randomization of photon positions. We remove 
this in our reduction of the data. In lStrickland et alJ (l2004l) we 
applied a technique called PHA randomization to remove an 
artificial quantization in the raw photon pulse height ampli- 
tude (PHA) values. CTI and time-dependent gain corrections 
automatically remove this quantization, making CIAO PHA 
randomization unnecessary (L. Townsley, private communi- 
cation). Removal of events at node-boundaries and standard 

Vikhlinin et al. , 2003, Corrections for time-dependence of ACI S gain, 
|http : //hea-www . harvard ■ edu/~alexey/acis/tgain/| 
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TABLE 1 

Chandra AND XMM-NEWTON OBSERVATIONS OF M82 



Date 


Observatoiy 


Instrument 


ObsID 


Exposure time 


Name used 


1999-09-20 


Chandra 


ACIS-1 


361 


15388 (13) 


(not used) 


1999-09-20 


Chandra 


ACIS-1 


1302 


32932 (13) 


33ks ACIS-I 


1999-09-20 


Chandra 


ACIS-1 


361-H1302 


48320 (13) 


merged ACIS-I 


2001-05-06 


XMM-Newton 


EPIC 


0112290201 


23327 (FN), 29685 (MOSl), 29687 (M0S2) 


XMM-short 


2002-06-18 


Chandra 


ACIS-S 


361 


17912 (S3) 


ACIS-S 


2004-04-21 


XMM-Newton 


EPIC 


0206080101 


59990 (FN), 73637 (MOSl), 73629 (M0S2) 


XMM-long 



Note. — The ObsID is the unique observation identification number assigned by the Chandra and XMM-Newton Science 
Centers. The exposure times (in seconds) quoted are those remaining after the removal of periods of higher-than-normal 
background. The specific chip (for Chandra observations) or detector (for XMM-Newton) the exposure time applies to is noted 
in parentheses. 



event grade filtering was performed to further reduce the num- 
ber of non-X-ray events in the data. 

Following gain or CTI correction, both the ACIS-S and 
ACIS-I datasets were processed to remove artificial streaks 
of events along the same CCD columns as the brightest point 
source in M82. This streak is caused by photons from bright 
point sources that are detected during the 40/is time period 
while the chip is being read out after each 3.2 s exposure. 
Corrections were applied to each ACIS data set to remove the 
event streaks expected from the six brightest point sources 
within the nuclear region of M82. 

In order to achieve the most-accurate absolute astrometric 
solution for these observations the latest satellite aspect so- 
lutions were applied to the data^. With these sub-arcsecond 
corrections applied, the astrometry of the data should be ac- 
curate to 0"6 at 90% confidence. 

2. 1.2. Flare-removal and source searching 

We used the techniques for removing periods of higher- 
than-normal particle background (flares), and point sourc e de- 
tection and removal, described in lStrickland et al.l (l2004l) . We 
refer the reader to that work for a detailed description. 

After flare removal, the exposure time on the S3 chip in the 
ACIS-S observation of M82 is te^p = 17.912 ks, texp = 32.932 
ks on the 13 chip of the 33ks ACIS-I observation and f„p = 
48.320 ks on the 13 chip of the merged ACIS-I observation. 

2.1.3. Background subtraction 

The diffuse hard X-ray emission from the central region of 
M82 covers a region of sky < 1' in angular diameter In this 
aperture the background (a sum of the extragalactic and galac- 
tic X-ray background emission and particle events in the de- 
tector) is weak but non-negligible. 

Soft diffuse X-ray emission from the superwind covers the 
major ity of the S3 chip i n the ACIS-S observation (see Fig. 3 
in Stri ckland et al.ll2004l) . leaving little area (< 5.4 arcmin^) 
from which to extract a high signal to noise background spec- 
trum. We instead used the blank-sky datasets produced by 
Maxim Markevitch^ for background subtraction. These in- 
clude CTI and gain correction consistent with our ACIS-S ob- 
servation of M82, unlike the blank sky datasets that are part 
ofCALDB 3.0. 

We used data from regions on the other ACIS chips that are 
free of point source or diffuse X-ray emission in our obser- 



' See |http : //asc .harvard, edu/ciao/ threads /arose c_cor 
* Markevitch. 2001, General discussion of the qui- 

escent and flare components of the ACIS background. 
http://cxc. harvard . edu/ contrib/maxim/bg/ index . html 
We used the 2004 December 12 version of the ACIS background datasets. 



vation (specifically chips SI, S2 and S4) to (a) check that the 
spectral shape of the background in our observation is consis- 
tent with that in the blank-sky datasets, and (b) calculate the 
factor by which the blank sky background data needs to be 
scaled to match the background in our observation. The back- 
illuminated S 1 chip is the most useful chip in this regard. In 
the ACIS-S observation the SI chip is offset to the SSE of 
the nucleus of M82, along the minor axis of the galaxy. The 
edge of the chip nearest to M82 is ~ 10' from the nucleus, 
well beyond the edge of any diffuse X-ray or other emission 
associated with the superwind. Furthermore, its sensitivity to 
background events is very similar to that of the S3 chip, unlike 
the front-ifluminated S2, S4 and S5 chips of the ACIS-S. We 
found that both the absolute normalization and spectral shape 
of the observed S 1 background matched that in the blank-sky 
background datasets very closely. The background normal- 
ization in our M82 observation was marginally lower than in 
the blank-sky datasets (by 1 .9 ± 1 .4% in the 5-8 ke V energy 
band, 1 .2 ± 1 .2% in the 2-8 ke V band, using an area of ^ 60 
arcmin^. These values are consistent the ~ 2% rms scatter in 
the normalization o f the hard X-ray component f rom observa- 
tion to observation (i Hickox & Markevitchll2006l) . Their mea- 
surements used a higher energy range than the one we have 
used, but as they also show that the spectrum of the back- 
ground very stable we conclude that a 1.9% renormalization 
is within the expected levels of variation. 

At lower energies, or in narrower energy bands the con- 
straints are weaker and hence consistent with no difference 
in normalization). There is known soft diffuse X-ray emis- 
sion at the location of the S2 and S4 chips, but in the hard 
X-ray band the observed and blank-sky background surface 
brightness levels in these chips are also consistent with each 
other. The resulting blank-sky background dataset for the S3 
chip had an effective exposure of 458.60 ks, when the 1.9% 
renormalization described above was been applied. 

The resulting background dataset has the same pointing and 
orientation as the actual observations of M82. Background 
subtraction is accomplished simply by subtracting from the 
observation data (be it an image or spectrum) the equivalent 
data from the background dataset in the same spatial region 
and energy band, scaled by the ratio of exposure times. The 
background is typically small compared to the true emission 
from M82 (see Table |2]i. At energies between E = 5 -1 keV 
the background within a radius of ^ 0.'5 of the center of M82 
is approximately 1 % of the total emission and ^ 4% of the 
-diffuse emission. 



^°iiiis method ignores the fact that some of the extragalactic 
X-ray background will be absorbed by M82's disk (i.e. M82 
will cast a shadow in the X-ray background), resulting in a 
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slight over-subtraction of the background. iBregman & IrwinI 
(|2002) demonstrated such shadowing of the soft X-ray back- 
ground by the disk of the edge-on spiral NGC 891. Shadow- 
ing is only visible at low energies (E <l keV), where the total 
Chandra ACIS-S background (instrumental particle events as 
well as the Galactic and extragalactic X-ray backgrounds) is 
reduced to ~ 70± 10% of the value found in regions not cov- 
ered by the galaxy. Thus we may have over-estimated the 
background in some regions by ^ 30%. Given that the back- 
ground is so small compared to the diffuse emission we are 
interested in this systematic error is negligible compared to 
statistical uncertainties in the diffuse flux (Table |2]i. Further- 
more the shadowing of the extragalactic background will be 
greatest toward the center of M82, but the relative significance 
of this systematic over-subtraction is mitigated by the intrin- 
sically higher emission from M82 in this region. 

For the ACIS-I observations of M82 background subtrac- 
tion is slightly more straightforward. Moderately large frac- 
tions of each of the four ACIS-I chips lie beyond the diffuse 
soft X-ray emission associated with the superwind, and hence 
are used to for background subtraction. All of the hard dif- 
fuse X-ray emission lies within the 13 chip, for which a back- 
ground region of area 18 square arcminutes free of diffuse and 
point-source X-ray emission is available. 

This local background estimate needs to be scaled appropri- 
ately when calculating the expected background in any other 
region on the same ACIS-I chip, given that the background 
will change with position. For image-based analysis we used 
the relative intensity from an exposure map (the product of the 
effective exposure, and collecting area, at some X-ray photon 
energy for each point on each ACIS chip) to map from an 
observed local background to a chip- wide model background 
image. This was done on a chip-by-chip and energy band-by- 
energy band basis. 

For spectral analysis we do not use this exposure-map based 
correction. We instead scale the spectrum from the back- 
ground region by the ratio of the area of the source region (the 
center of M82) to the larger background region. No correction 
for vignetting or shadowing of the background is applied, but 
this does not introduce significant systematic errors because 
the background is so small compared to the diffuse emission. 

2.1.4. Spectral fitting 

Spectral fitting was performed using XSPEC (version 
11.3. It). Spectral response files were created for the ACIS-S 
data using CIAO and the latest CALDB. The instrument effec- 
tive area takes into account the time-and-spatially-dependent 
decrease in low-energy sensitivity due to molecular contami- 
nation of the optical blocking filter (OBF), as well as the weak 
spatial dependence of the effective area. For the ACIS-I ob- 
servation we used the ACIS -13 aim-point RMF file provided 
as part of the PST CTI correction software, and created the 
ARF file with the CIAO task MKARF using the QEU file pro- 
vided with the PSU software. 

For the ACIS-S observation a background spectrum was 
extracted from the blank-sky background dataset in a 1'- 
diameter circle centered on the dynamical center of M82 (see 
§|4]i. Note that using the blank-sky background datasets auto- 
matically takes account of the spatial dependence of the back- 
ground. For the ACIS-I observation the background spectrum 
was taken from the large emission-free region on the 13 chip 
mentioned above. 

Background-subtracted spectra were manually created us- 
ing the CIAO task DMTCALC, instead of letting the XsPEC 



spectral-fitting program perform the background subtraction 
itself. This allows us to both re-bin the spectra to ensure a 
minimum of 10 counts per bin after background subtraction 
(using DMGROUP), an d associate mo re-realistic statistical er- 
rors with the data (the lGehreli (Il986l) approximation to Pois- 
sonian errors, instead of ^/N errors). 

2.1.5. The ejfect of new calibration products 

We investigated the effect these recent calibrations have 
on our results, compared to those based on the CIAO 2.2, 
CALDB 2.12, non-C TI-corrected data reduction used in 
[Strickland et al.l (l2004l) . There is no significant difference in 
conclusions based on image analysis using these different cal- 
ibrations. This is to be expected, as the photon energy range 
in most images is much larger than the changes to photon en- 
ergy resulting from the CTI or gain corrections. Where the 
new calibrations do make a significant difference is in spec- 
tral fitting, in particular to the energy of line features or line- 
complexes. For example, we find the mean energy of the 
Fe-K feature in the ACIS-S observation nuclear spectrum is 
~ 0.1 keV higher when using the most recent calibration data 
and gain corrections. 

2.2. XMM-Newton observations 

The two longest XMM-Newton observations of M82 were 
performed on 2001 May 06 (ObsID 112290201, PI Mar- 
tin Turner, initial duration 29.4 ks, henceforth referred to as 
"XMM-short") and on 2004 April 21 (ObsID 206080101, ini- 
tial duration 105.4ks, PI Piero Ranalli, which we shall refer to 
as "XMM-long"). We did not use a shorter observation (^ 10 
ks) taken on the same date as the XMM-short observation be- 
cause of the low S/N of the resulting spectra. 

We reprocessed the archived data using version 6.5 of the 
XMM-Newton SAS software, and reran the EPIC MOS and 
PN analysis chains with the latest calibration data. The 
recommended standard event grade filters and proton flare 
screening was applied to each dataset. After filtering the re- 
maining exposure times in the XMM-short observation were 
2 1 . 1 ks and 29.4 ks for PN and each MOS instruments respec- 
tively, and in the XMM-long observation 54.0 ks and 72.8 ks 
for the PN and each MOS detectors respectively. 

Spectra and response files were created using SAS us- 
ing the methods described in the on-line SAS Data Analysis 
Threads^. 

Nuclear-region spectra were extracted from the same area 
used for the Chandra ACIS spectra (see § 14. 11 1. Background 
spectra for each observation were accumulated in regions of 
radius 1 .'25 that were free of point sources or diffuse X-ray 
emission. The same background region was used for all three 
detectors (PN, MOS 1 and MOS2) and whenever possible was 
on the same CCD chip as the nuclear diffuse emission. 

3. IMAGE ANALYSIS OF THE CHANDRA DATA 

3.1. The spatial distribution of the diffuse soft and hard 
X-ray emission 

The spatial distribution of the soft and hard X-ray emission 
within ^ 1 kpc of the center of M82 as seen in the ACIS-S 
and merged ACIS-I observations is shown in Fig. [T] Point- 
like sources identified within each observation are enclosed 
by circular or rectangular apertures. The only significant dif- 
ferences between the ACIS-S and ACIS-I images are changes 

' See |http : //xmin. vilspa ■ esa ■ es/sas/new/documentation/threads/l 
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Fig. 1. — The center of M82 as seen with Chandra. Each image shows the raw, unsmoothed, photon distribution in the soft 0.3-2.0 keV (panels a and c) 
and and hard 2.0-8.0 keV energy bands (panels b and d), binned in 0"492-wide pixels, over a 2 X 2 kpc (lf90 X lf90) region centered on the dynamical center 
of the galaxy. Panels a and b are from the ACIS-S observation, while panels c and d are from the merged ACIS-I observation. The images are shown on a 
square-root intensity scale running from to 75 counts per pixel. Point-like X-ray sources detected in the soft and hard energy bands are shown surrounded by 
a circle, equivalent in size to the region used to remove the point sources from the images and spectra. Panels e, f, g and h are equivalent to panels a, b, c and d 
respectively, except that point sources have been removed and i nterpol ated over, and the intensity scale runs from to 25 counts per pixel. The background has 
been subtracted from all the images in the manner described in § 12.1.31 



in the brightness of some of the point-like X-ray sources due 
to intrinsic variability, and drops in the apparent diffuse soft 
X-ray surface brightness in the bottom and top left of the 
ACIS-I image due to the gaps between the CCD chips. 

Images with the point-sources removed are shown in panels 
e to f of Fig. [T] and were created by replacing events within 
the point source regions by a Poissonian random deviate con- 
sistent with the events detected within a local background an- 
nulus of width l."5 around each source-region (CIAO task 
DMFILTH run in POISSON mode). Care was taken to ensure 
that these background regions did not themselves include any 
other point source emission. This source removal method is 
performed on each image. Images in an energy band < 3 keV 
have had the point sources removed that were detected only 
in the soft E = 0.3-2.0 keV energy band, as sources detected 
only within the hard E = 2.0-8.0 keV energy band are most 
likely heavily obscured and hence invisible at lower energies. 
Hard X-ray images are cleared of point sources detected in 
either of the soft or hard X-ray bands. This method produces 
a source-subtracted image of the 33ks ACIS-I data that is es- 
sentially identical that from the PSF-fitting and subtraction 
method presented in Griffiths et al. (2000). 

A more detailed break down of the spatial distribution of 
the diffuse X-ray emission within the central 2 kpc of M82 
as a function of photon energy is shown in Fig. |2] based on 
the ACIS-S data alone. For comparison we also plot ground- 
based optical R-band and continuum-subtracted Ha-i-[N II] 
images of the same region. 

The effect of absorption by material between us and in the 
X-ray emitting plasma near the plane of the galaxy is visible at 
X-ray energies £ < 1 . 1 ke V as a band of low X-ray surface- 
brightness bisecting the wind from north east to south west, 
obscuring the nuclear starburst region. This X-ray absorption 
band has a larger minor-axis extent than the narrower dust 



lanes seen in the optical images shown in Fig. |2] 

At X-ray energies E < 0.6 keV (at energies around the 
E=0.57 O VII and 0.65 keV O VIII emission Hnes) the cen- 
tral region of the galaxy is heavily obscured, and the bright- 
ness of the diffuse X-ray emission peaks at z > 20" (~ 400 
pc) from the mid-plane. In contrast the mid-plane, nucleus 
and the inner wind are visible in emission from Helium 
and Hydrogen-like ions of Magnesium (E ^ 1.3 - 1.5 keV, 
where the similarity in morphology to the Ha-i-[N II] emis- 
sion is striking) and Silicon (£ ~ 1.8 to 2.0 keV), and the 
brightest emission occurs within z ~ 50 pc of the mid-plane. 
Read & Stevens ( 2002) and Origlia et al. (2004) studied the 
nuclear-region element abundance (within ^15" radius of the 
cent er of the gala xy) using XMM-Newton RGS grating data. 
The lOrigha et al.l study derived a hot-phase Oxygen abun- 
dance that is unusually low compared to other a-elements. 
Neither study could take account of the strong spatial varia- 
tions in foreground absorption within this region (because of 
the low spatial resolution of XMM-Newton), which our Chan- 
dra data shows to be of particular significance in M82 at en- 
ergies near the oxygen lines. Fitting simple models to spec- 
tra with intrinsically complex absorption is known to lead to 
systematic errors in deriv ing X-ray abundanc es from CCD- 
resolution X-ray spectra (IWeaver et al.l 120001) . We therefor 
suspect that the abu ndances derived from X MM- Newton stud- 
ies of M82 in the iRead & StevensI (12002*) and 'O riglia et all 
(2004), in particular the apparently low Oxygen abundance, 
are incorrect. 

Diffuse hard X-ray emission is clearly visible within the 
central starburst region, even in images without point source 
subtraction. Unlike the diffuse soft X-ray emission from the 
superwind, the spatial extent of the brightest region of diffuse 
hard X-ray emission is larger along major axis (i.e. along the 
plane of the galaxy) than along the minor axis. Measured 
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by eye, hard diffuse X-ray emission has a major-axis radial 
extent of ~ 20" 350 pc) approximately centered on the 
dynamical center of the galaxy, and extends ^ 8" (~ 140 pc) 
and ^ 12" (^^ 210 pc) above and below the plane of the galaxy 
(to the south east and north west of the plane, respectively). 
Note that diffuse hard X-ray emission can be found outside 
this central region, although at lower surface brightness. This 
is discussed further in § 13.31 

The spatial distribution of the diffuse hard X-ray emission 
is relatively smooth. This is not a limitation of spatial resolu- 
tion, photon statistics or the method of point source removal, 
as significant structure is visible in the soft X-ray images. This 
suggests that the source of the hard X-ray emitting material is 
also relatively uniformly distributed within the starburst re- 
gion. 

There are three bright point-like X-ray sources just to the 
west of the center of the galaxy in the ACIS-S and ACIS-I 
observations (one of wh ich is the variable ultra-luminous X- 
ray so urce described in iKaaret et al.ll200Tt iMatsumoto et alJ 
1200 Ih . Point-source removal is difficult in this region, as 
bright sources have large PSF wings. Combined with the 
close proximity of these sources to each other we were forced 
to use a 10."8 x 4."4 rectangular region for point source re- 
moval. The flux at each pixel within this region in the source- 
subtracted images is a Poisson random deviate drawn from the 
distribution of the fluxes within 1 ."5 of the edge of this rectan- 
gle, which includes the region of brightest diffuse hard X-ray 
emission at the mid-plane of the disk. This method is likely to 
slightly over-estimate the diffuse X-ray surface brightness at 
the off-plane location of the ULX, as this entire region is filled 
in by DMFILTH approximately uniformly. Thus the apparent 
association between the ULX and a bright region of diffuse 
hard X-ray emission extending off the bright mid-plane ridge 
(e.g. seen in Fig.|2ji) may not be real. 

3.2. Diffuse emission and unresolved point sources 

Count rates for the diffuse emission and the detected point 
sources within 500 pc of the center of M82 are given in Ta- 
ble |2] for both the ACIS-S and merged ACIS-I observations 
(the results for the 33 ks ACIS-I observation are consistent 
with the merged observations). The diffuse count rates are 
based on the point-source removal technique described above. 
The point source count rate within any energy band is the dif- 
ference between the total count rate and the estimated diffuse 
count rate. 

The fraction of the X-ray counts within each energy band 
from apparently diffuse emission, the diffuse fraction /diff, is 
also presented in Table|2] Some fraction of the apparently dif- 
fuse emission comes from unresolved point sources, as only 
the members of the M82 X-ray source population with in- 
dividual luminosities Lx > 4 x 10^*^ erg s"' are detected as 
point-like in these observations dStrickland et al.l i2004). Star 
forming galaxies typically have X-ray point source popula- 
tions with cumulative lu minosity functions of the form A^(> 
L) oc L"', where 7 < 1 ( Colbert et an |2004)). i.e. the total lu- 
minosity is dominated by the brightest sources. Thus the ratio 
of the total luminosity in detected point sources to the total 
luminosity from all point sources is x,- w 1 - (imin/i'max)' 
Here L^in and Lmax are the luminosities of the faintest and 
brightest point sources detected. Using this to correct the dif- 
fuse fraction for the expected unresolved point source con- 
tribution we obtain an estimate of the true diffuse emission 
fraction /Jjjj- = 1 - (1 -fduij/xr. We used a variety of methods 
to estimate Xr in the ACIS-S and ACIS-I observations of M82, 



where the faintest point sources detected in any energy band 
typically have > 10 counts, and the brightest source is the 
ULX with several thousand counts. We found that for M82 Xr 
depends only weakly on the energy band considered, and the 
largest effect is the assumed luminosity fun ction slope. We 
used two val ues for the slope 7 = 0.50 (Kilgard et al.ll2003h 
and 7 = 0.57 dColbert et alj|2004|) . deriving estimates of Xr in 
the range 0.87 - 0.93 for the ACIS-S observation, and Xr in the 
range 0.91 - 0.96 for the deeper merged ACIS-I observation. 

Correcting /diff for unresolved point sources has little effect 
at energies below E ^ 3 keV, but is a non-negligible correc- 
tion at higher energies. Nevertheless, significant residual dif- 
fuse hard X-ray emission exists even after correcting for un- 
resolved point sources, with a fraction > 70 % (/^ff//diff) of 
the apparent diffuse hard X-ray emission being unattributable 
to unresolved point sources. 

In general the estimated diffuse fraction is consistent be- 
tween the ACIS-S and ACIS-I observations, despite the dif- 
ferences in spectral sensitivity between the back-illuminated 
and front illuminated chips, and the different epochs of obser- 
vation. Note that the diffuse fraction is likely to vary with 
time, in particular in the harder X-ray bands {E > 3 keV) 
as significant variabiUty is associate d with the nuclear X- 
ray poi nt sources (e.g. [Co llura et anil994t [P tak & Griffith^ 
1999; Matsumoto & Tsurul ll999t ~ lMats'umoto et al.i i20o!t 
Kaaret et al. 2001). Several of the point-like X-ray source did 
show minor variability in flux between the ACIS-I and ACIS- 
S observations (as can be seen in Fig.[T]|, but this was not of a 
magnitude sufficient to significantly alter the diffuse fraction. 

The diffuse fraction is high (/Jff > 70%) at energies E <3 
keV, but is relatively constant with energy at higher energies 
(£■ > 3 keV) at a value of ^ 20 - 30%. This transition in 
diffuse fraction occurs at the same energy associated with the 
transition from X-ray emission extended preferentially along 
the minor-axis (the soft X-ray emission in the wind) to pre- 
dominantly major-axis extended hard diffuse X-ray emission 
shown in Fig.|2] 

Although the ULX is the single brightest X-ray point source 
over the 0.3-8.0 keV energy band, it accounts for somewhat 
less than half of the total resolved point source related flux 
from the central kiloparsec of M82 in either Chandra ACIS 
observation. In the merged ACIS-I observations of 1999 
September 20, the ULX contributed - 51% of the 0.3-2.8 
keV energy band count rate, ~ 41% of the 3.0-9.9 keV band 
count rate and ^ 45% of the 0.3-9.9 keV band count rate 
within a radius of 29" of the center of the galaxy. In the ACIS- 
S observation of 2002 June 18 the corresponding numbers are 
19, 39 and 30% respectively. We estimate that the true count 
rate of the ULX is probably 50% higher than the observed 
count rate in the soft energy band, and 20% higher in the hard 
energy band, due to the effects of pile-up in the ACIS CCD 
detectors (based on the PIMMS count rate prediction tool), 
but this does not alter the general conclusion that in total the 
other resolved point sources contribute a comparable flux to 
that of the ULX at the specified epochs. 

3.3. The spatial extent of the diffuse hard X-ray emission 

We created surface brightness profiles along the major and 
minor axes of the galaxy in several different energy bands in 
order to better quantify the spatial extent of the diffuse hard 
X-ray emission in M82, particularly in comparison to the bet- 
ter studied diffuse soft X-ray emission. We also calculated 
the spectral hardness Q = (H-S)/(H+S), where H and S are 
the counts in chosen hard and soft energy bands respectively. 
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TABLE 2 

Count rates and diffuse emission fractions. 



Energy 
(keV) 
(1) 



^tot 
(cts/s) 
(2) 



(cts/s) 

(3) 



(cts/s) 
(4) 



Rhg 

(cts/s) 
(5) 



/diff 

(6) 



-'dift 
(7) 



ACIS-S observation (S3 chip) 



0.3 
0.6 
1.1 
1.6 
2.2 
0.3 
5.0 
6.0 
3.0 



-0.6 

- 1.1 

- 1.6 
-2.2 
-2.8 
-2.8 
-6.0 
-7.0 
-9.9 



0.0208 ±0.0012 
0.5603 ±0.0057 
0.4783 ±0.0052 
0.3298 ± 0.0044 
0.1256 ±0.0027 
1.5 149 ±0.0093 
0.0485 ±0.0017 
0.0289 ±0.0013 
0.3166 ±0.0043 



0.0191 ± 
0.5263 ± 
0.4347 ± 
0.2727 ± 
0.0901 ± 
1.3429 ± 
0.0138 ± 
0.0079 ± 
0.1071 ± 



0.0011 
0.0055 
0.0050 
0.0040 
0.0023 
0.0087 
0.0009 
0.0007 
0.0026 



0.0017 ±0.0004 
0.0340 ±0.0014 
0.0437 ±0.0016 
0.0571 ±0.0018 
0.0356 ±0.0015 
0.1721 ±0.0032 
0.0347 ±0.0014 
0.0210 ±0.0011 
0.2095 ±0.0035 



0.00089 ± 0, 
0.00051 ±0, 
0.00030 ±0, 
0.00057 ± 0, 
0.00025 ± 0, 
0.00252 ±0, 
0.00037 ±0, 
0.00041 ±0, 
0.00548 ±0, 



.00005 
.00004 
00003 
00004 
.00003 
.00008 
00003 
.00003 
.00011 



0.917 ±0.203 
0.939 ±0.041 
0.909 ±0.035 
0.827 ±0.029 
0.717 ±0.033 
0.886 ±0.017 
0.284 ±0.015 
0.274 ± 0.020 
0.338 ± 0.007 



0.91 ±0.20 

0.93 ±0.05 
0.90 ±0.04 
0.81 ±0.04 
0.69 ± 0.04 
0.87 ±0.03 
0.20 ±0.01 
0.19 ±0.02 
0.26 ±0.01 



Merged ACIS-I data (13 chip) 



0.3- 
0.6- 
1.1- 
1.6- 
2.2- 
0.3- 
5.0- 
6.0- 
3.0- 



0.6 
1.1 
1.6 
2.2 
2.8 
2.8 
6.0 
7.0 
9.9 



0.0085 ± 0.0005 
0.3252 ±0.0026 
0.3773 ±0.0028 
0.2276 ± 0.0022 
0.0949 ±0.0014 
1 .0336 ± 0.0047 
0.0405 ± 0.0009 
0.0233 ± 0.0007 
0.2473 ± 0.0023 



0.0082 ± 
0.3172± 
0.3447 ± 
0.1877± 
0.0691 ± 
0.9269 ± 
0.0101 ± 
0.0063 ± 
0.0830 ± 



0.0005 
0.0026 
0.0027 
0.0020 
0.0012 
0.0044 
0.0005 
0.0004 
0.0014 



0.0004 ± 
0.0080 ± 
0.0326 ± 
0.0398 ± 
0.0258 ± 
0.1067± 
0.0305 ± 
0.0170 ± 
0.1644± 



0.0001 
0.0004 
0.0008 
0.0009 
0.0008 
0.0015 
0.0008 
0.0006 
0.0019 



0.00035 ±0, 
0.00038 ± 0, 
0.00022 ± 0, 
0.00025 ± 0, 
0.00020 ± 0, 
0.00141 ±0, 
0.00024 ± 0, 
0.00024 ± 0, 
0.00212 ±0, 



00011 
.00011 
.00009 
00010 
.00009 
.00019 
00009 
00009 
.00023 



0.956 ±0.288 
0.975 ± 0.053 
0.914 ±0.025 
0.825 ±0.021 
0.728 ±0.024 
0.897 ±0.013 
0.249 ± 0.009 
0.269 ±0.013 
0.335 ±0.005 



0.95 ±0.23 
0.97 ±0.06 
0.91 ±0.04 
0.81 ±0.03 
0.71 ±0.03 
0.89 ±0.03 
0.20 ±0.01 
0.22 ±0.01 
0.29 ±0.01 



Note. — Column 1 : Energy band associated with the following count rates. Columns 2, 3 and 4: Background-subtracted 
count rate within a 0f476 radius of 09:55:51.9, +69:40:47.1 (J2000.0). The total background-subtracted count rate (including 
both diffuse emission and point sources) is Stot, while Rjiff and i?ps are the count rates associated with d iffuse and point- 
source emission respectively. Column 5: The background count rate within the same region (see § 12.1.3) . Column 6: The 
apparent diffuse fraction f^nf = Rduf/Rtat- Column 7. The diffuse fraction after correcting for unr esolv ed point sources, 
assuming x,- = 0.90 for the ACIS-S observation and jr,- = 0.94 for the merged ACIS-I observation. See j} l3.2l for further details. 
The errors quoted in columns 2-6 are 1 a statistical uncertainties, while the error quoted in column 7 is the combination of 
the statistical uncertainty in /diff and an assumed 3% uncertainty in .v, . 



along these axes. Only data within 48" (835 pc) of the re- 
spective axis was used when creating each profile in order to 
achieve a reasonable balance between maximizing the num- 
ber of counts without averaging over too large an area. We 
re-binned each profile so as to achieve > 20 counts per bin 
per energy band for plotting and fitting the surface brightness 
profiles, and > 40 counts per bin per energy band when cal- 
culating the spectral hardness. The resulting profiles for the 
merged ACIS-I observation is are shown in Fig. [3] The ACIS- 
S surface brightness profiles are in general very similar to the 
ACIS-I profiles shown in this figure. The most notable dif- 
ferences are purely instrumental due to the effects of the gaps 
between the ACIS-I CCD chips (at z -600 pc and r - -400 
pc). At photon energies E > I keV all profiles show a cen- 
tral peak in emission, the surface brightness of which drops 
off rapidly with distance, surrounded by more-extended lower 
surface brightness X-ray emission. 

Along both major and minor axes the hardness ratios show 
that the emission within < 400 pc of the center of the galaxy 
is spectrally distinct from the kiloparsec scale emission. The 
effect of absorption by the disk is seen most clearly in the mi- 
nor axis hardness ratios. Along the minor axis this hardness 
ratio drops from the nuclear starburst region to the larger scale 
wind. Moving to the SE along the wind the hardness is rela- 
tively constant in both the ACIS-S and ACIS-I data, although 
in contrast the hardness ratio increases moving NW from the 
nucleus. 

Simple model fitting demonstrates that to first order the 
scale heights of the soft and hard diffuse X-ray surface bright- 
ness outside the nuclear region are very similar, H ^ 500 - 
1000 pc. A particularly interesting result is the apparent de- 



tection of diffuse hard X-ray e mission out to l arger z and r 
distances than recognized in the [Griffiths et alJ (l2000h analy- 
sis of the initial ACIS-I observations (out to \z\ ^ |r| 2 kpc). 
However this faint extended hard emission X-ray might be 
due to the extended wings of the Chandra point spread func- 
tion, rather than being genuinely diffuse emission. We will 
investigate this issue further in the future. 

3.4. Comparison to features at other wavelengths 

It is instructive to compare the spatial distribution of 
the soft and hard X-ray within the central region of 
M82 to starburst-related activity seen at other wave- 
lengths. In Fig. |4] we plot the locat ion of the near-IR 
super-star c lusters (McCrady et al. 2003), the radio SNRs 
(lUnger et a l., 1984: Kronberg & Sramek 1985; Muxlo w et all 
119941: IWills etal.1 1 1998b. and the ex panding shells seen 



descri bed in Hi dWills et alJ [2q02h and CO-emission^ 
(Neini nger etaP Il998t IWeiBet alJ Il999t [Matsushita et"aD 

2000), 

Matsumoto et al.l (1200 ll) discuss the relationship between 
the X-ray point sources seen in the Chandra High Resolu- 
tion Camera observations o f M82 and the 23 rad io sources 
seen at 5 GHz presented in iMuxlow et alJ (|1994|) . The co- 
ordinates of the X-ray point sources seen in both our re- 
processed ACIS observations and the published HRC obser- 
vations match within the uncertainties. In Fig. |4]we plot the 
locations of the 33 compact radio sources observed at 1.4 



* IWills et alj <2002l) show that their shell 3 is most-probably the 
northern side of the C O-superbubble discussed in WeiB et al. ( 1999) and 
[Matsushita et alj j2000h . the exact energetics and mass of which are not 
agreed upon. 
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Fig. 2. — The distribution of diffuse X-ray emission in the central 2x2 kpc of M82 as a function of energy is shown in panels a to g, compared to optical 
R-band continuum and continuum-subtracted Hq+[N II] emission in panels h and i. The energy bands used are E=0.3-0.6 keV (a), 0.6-1.1 keV (b), 1.1-1.6 keV 
(c), 1.6-2.2 keV (d), 2.2-2.8 keV (e), 3.0-7.0 keV (f) and 6.0-7.0 keV (g). All images are shown on a square root intensity scale. Contour levels increase in 
factors of 2 in surface brightness. Images are aligned so that North is up and East is to the left. The X-ray images are from the ACIS-S observation, have had 
background and point source emission removed, and have been adaptively smoothed so that features are locally significant at S/N=3. The intensity scale for the 
X-ray images runs from to 5 counts (per Of'492 pixel), and the contours begin at a level of 0.04 counts per pixel. The effect of strong absorption in the plane of 
the galaxy is visible at X-ray energies < 1 . 1 keV as a band bisecting the wind from north east to south west, obscuring the nuclear starburst region. 



GHz by IWills et aH (|1998h . The clearest association of an 
X-ray source with a radio source at either 1.4 or 5 GHz is 
with the enigmatic radio source 41.9H-58, on e of the few ra- 
dio so urces in M82 that is fading with time (iKronberg et alJ 
llOOOV The radio source 44.0H-59.6, another of the bright- 
est radio sources, falls at the location of a weaker hard X- 
ray point source. In general the radio sources, most of which 
are yo ung SNRs evolv ing in co nfining dense gaseous env iron- 
ments (IKronberg et a l. 2000; C hevalier & Franssonl2001h . are 
not associated with X-ray point sources. 

Although the SNRs associated with radio sources may not 
directly contribute to the mechanical powering of the super- 
wind because of their confinement, they must delineate the 



region associated with current SN activity. It is noteworthy 
that the major-axis extent of both the soft and hard diffuse X- 
ray emission is so similar to that of the radio sources — this 
is the unambiguous sign of a massive-star-related SN-driven 
wind. The distribution of the diffuse hard X-ray emission dif- 
fers from that of the radio sources only in that it has a larger 
minor-axis extent, to be expected if the X-ray-emitting plasma 

is being advected out along the minor axis. 

The near-IR star clusters defined in iMcCrady et alJ (l2003h 
are more widely distributed along the major axis of the galaxy 
than the radio sources or the soft and hard diffuse X-ray emis- 
sion. However not all of these clusters are young enough to 
support ongoing core-collapse SN activity. It is not surprising 
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Fig. 3. — Minor and major axis diffuse X-ray emission surface brightness and sp ectral h ardness profiles, from the merged ACIS-I observation. The positive z 
axis is along PA = 343°, and the positive r axis along PA = 253°. The background (§ l2.1.3t has been removed, and point sources interpolated over. Error bars are 
la uncertainties. The spectral hardness ratio is defined as Q = {H -S)/{H + S), where H and S are the counts in the higher and lower energy bands respectively. 

spatially-varying absorption and possible variations in mean 
plasma temperature occur (see Fig. [Sj. With the AE ^ O.l 
keV resolution of the AClS-1 and ACIS-S CCD detectors the 
strong line emission from O, Fe, Ne, Mg, Si, and S are heav- 
ily blended, making a unique and physically accurate spectral 
fit currently difficult. Our focus in this paper is instead on 
the harder X-ray emission with energies E >2 keV, where 
there are fewer spectral features and absorption can effec- 
tively be ignored. In particular we are interested in the con- 
tinuum spectral shape of the diffuse hard X-ray emission, and 
whether it is similar to that of the resolved point source pop- 
ulatio n. Given the am biguous detection of iron line emission 
in the [Griffiths et al.l ( 2000) study, another question we wish 
to answer is what are the statistical significances and fluxes 
of any E ^ 6.4, 6.7 or 6.9 keV iron line emission from the 
diffuse component? Is the Fe-K emission associated with the 
point source pop ulation, and/or the enigmatic lu minous and 
variable ULX, as IStrohmaver & Mushotzkvl (12003 ) conclude 
based on the ObsID 1 12290201 XMM-Newton observation of 
M82? 



that there is very little diffuse X-ray emission in t he vicinity of 
the 60 ± 20 Myr-old super star cluster F (Smith & Gallaghe^ 
I2001h . One of the point-like X -ray sources may be associated 
with the nearby cluster L. 

Three o f the four expanding Hi shells identified by 
IWills et all (l2002h he within the region covered by diffuse 
soft and hard X-ray emission. As mentioned previously, one 
of these Hi s hells is probably a ssoc iated with the CO -bubble 
discussed in IWeiB et al] (119991) and [Matsushita et alj (I200Q) . 
This entire region along the major axis is the location of the 
brightest soft and hard diffuse X-ray emission, there being no 
discernible difference in X-ray surface brightness between the 
regions interior and exterior to the Hi and/or CO shells. 

4. SPECTRAL ANALYSIS OF THE CHANDRA AND XMM-NEWTON 

DATA 

We will begin with an analysis of the Chandra ACIS-S 
and ACIS-I spectra of each of the diffuse, point source and 
ULX spectra, followed by assessment of the significance of 
the E ~ 6.4 keV and E ^ 6.7 keV iron line features in those 
spectra. We will then present a more elaborate model that at- 
tempts to account for the presence of unresolved point source 
emission in the diffuse spectrum, and assess what limits can 
be placed onE = 6.97 keV emission from hydrogen-like ion- 
ized iron in the Chandra spectra. This section is concluded 
with analysis of the integrated nuclear spectra as seen in the 
two separated XMM-Newton observations (by integrated we 
mean that the spectra include diffuse, point source and ULX 
emission). We will show that a full picture of the origin of the 
iron line emission from M82 is only obtained if the data from 
both Chandra and XMM-Newton at all epochs is assessed. 

The Chandra ACIS spectra of the central kiloparsec of M82 
are quite complex, particularly in the soft X-ray band, where 



4. 1 . Chandra ACIS spectra 

We extracted spectra from the ACIS-S and ACIS-I datasets 
within a radius of 28 "6 (equivalent to a physical radius of 500 
pc) of the dynamical center of M82 at a = 09''55"51?9,(5 = 
+69°40'47"1 (J2000.0). We extracted and fit spectra to both 
the merged ACIS-I observation (which has the greatest num- 
ber of counts) and the 33ks ACIS- I observation (the obser- 
vation used by [Griffiths et alj[2000h . A spectrum of the dif- 
fuse X-ray emission was created by excluding all detected 
point sources using either circular or rectangular masks, as 
shown in Fig.[T] The source extraction radius used depended 
on source flux, and ranged from 1"2 to 3."0, most typically 
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Fig. 4. — The X-ray emission from the central 1.5 X 1 .5 kpc region of M82 is shown in c omparison to the n ear-IR super-star clusters of lMcCradv et al.l 120031 
small blue crosses), the radio SNRs of'Wills et al.' ('1998, red X's), the Hi shells described in lWills et alj j2002l shown as blue circles of diameter equivalent to the 
observed shell) and the CO-superbubble I WeiB et al. 1999, cyan polygon). Panels a and b show ACIS-S images in the E = 0.3-2.0 keV energy band, including 
the point source emission in (a) and source-subtracted in (b). Adaptive smoothing has been used in (a), whereas in (b) the smoothing used a Gaussian mask of 
FWHM= 2". Panel (c) is equivalent to (a), but shows the merged ACIS-I image. Panels (d) through (f) are equivalent to panels (a) through (c), except the images 
are in the 2.0 - 8.0 keV energy band. 



being 1 "5 - 2". These radii were chosen to maximize the re- 
moval of point-source related flux from the diffuse spectrum: 
a radius of 1 ."5 encloses 85 - 97% of the flux from a point 
source, depending on the energy of the photons (the enclosed 
energy fraction is smaller at higher energies). We also created 
a summed spectrum of the point source population (excluding 
the ULX) using these point source extraction regions, and a 
spectrum of the ULX itself. A small component of the point 
source and ULX spectra is from truly diffuse emission within 
the extraction regions. The spectrum of the ULX suffers from 
"pile-up," in which two or more X-ray photons arrive in a sin- 
gle pixel within a single 3.2s ACIS exposure and are falsely 
detected as a single photon of higher energy. We do not seek 
to correct for this process, but merely to see if any Fe-K line 
emission can be detected above the local continuum. 

We only fit the ACIS spectra in the energy range £ = 3.15 
-3.80 keV and 4.20 - 9.0 keV. Excluding the data below 3.15 
keV removes any contribution from the soft diffuse emission 
and associated Si and S emission, as well as the instrumental 
Ir edges. We exclude the energy range 3.8-4.2 keV to avoid 
any He-like and H-like Calcium emission, which is present 
at marginal significance in the ACIS-S diffuse emission spec- 
trum but is not obvious in the AClS-1 spectra. Our aim is to 
simplify the spectral fitting by essentially only considering a 
pure continuum plus any iron lines. 

When modeling the continuum in the diffuse and point 
source spectra we found that standard thermal bremsstrahlung 
and power law models were adequate. The spectra of the ULX 



are heavily piled up, and can not be adequately fit with power 
law or bremsstrahlung models (even when we also included 
the "pileup" model available in XSPEC). We found that a sim- 
ple broken power law model, with the transition between the 
two power law slopes occurring at E ^ 6.0 ±0.5 keV, pro- 
vided the best fit the the apparent continuum in the ULX spec- 
tra. No absorption model was used, as a column A^h ~ 10 is 
required to produce an optical depth of r '--^ 1 at £ ^ 3 keV. 

At the spectral resolution of ACIS-S and AClS-1 the Fe-Ka 
fluorescence lines from nearly neutral iron and the Fe Hea 
line complex from highly ionized iron can be modeled as nar- 
row Gaussian (cr = 1 eV) lines at E = 6.40 and E =6.69 keV 
respectively. Although it is possible to fit the ACIS spectra 
to determine the the line energy and/or width this does not 
place any useful constraint on these parameters, nor does this 
result in statistically superior fits compared to models using 
the expected line width and energy (by statistically superior 
we mean a change in chi-squared per change in degrees of 
freedom Ax^/Az^ > 1). For example, when fitting for the 
line energy of the Fe Hea line the best fits were E = 6.66^ \^ 
keV (ACIS-S) and £ = 6.69 ± 0. 17 keV (merged ACIS-I spec- 
trum). Thus for the majority of the spectral fits described be- 
low, the line energies and widths were fixed at the expected 
values, and we fit only for each line intensity. 

The best-fit spectral models are plotted along with the dif- 
fuse, point source and ULX spectra in Fig.|6] and the best-fit 
model parameters are shown in Table [3] The derived fluxes 
for the diffuse emission are shown in Table |4] The pure con- 
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TABLE 3 

Model fits to the diffuse, point source and ULX spectra 



Model Dataset Continuum model 6.4 keV Fe Ka line 6.7 keV Fe Heo line 6.9 keV Fe Lya line ^ 

kT or r norm norm EW norm EW norm EW 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 



Chandra diffuse emission spectra 



BR, two lines 



BR, no lines 



PL, two lines 



PL, no lines 



BR+PL, three lines 



PL+PL, three lines 



ACIS-S 
33ks ACIS-I 
merged ACIS-I 
ACIS-S 
33ks ACIS-I 
merged ACIS-I 
ACIS-S 
33ks ACIS-I 
merged ACIS-I 
ACIS-S 
33ks ACIS-I 
merged ACIS-I 
ACIS-S 
33ks ACIS-I 
merged ACIS-I 
ACIS-S 
33ks ACIS-I 
merged ACIS-I 



6 2+^-' 
"■•^-2.0 

3 4+0.7 

8.0:|| 
44+1.3 
4 

Z.1_Q3 

^■'-0.3 
2 7+"-^ 

2 0+^'^ 
^-^-0.2 

3.8+2-5 

2.8!*| 

7 7+0.6 
^■'-0.4 
2 5+"^ 

M 

] 

-0.2 



3.0; 

3.0 



9 7+2.8 
^•'-1.9 
13.0!35 

13.6!3:« 

8.8!2| 
10.9+f^ 
11.3+^' 
13.7«' 
26.4;'- 
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14.9L3 
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65.1 
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83 
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115 




87.2 


136 




41.9 


83 




71.4 


117 




105.4 


138 


37 


35.1 


80 


20 


63.1 


114 





91.0 


135 


24 


34.3 
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59.9 


114 





85.1 
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Chandra summed point source spectra, excluding the ULX 



, two lines 



, no lines 



ACIS-S 
33ks ACIS-I 
merged ACIS-I 
ACIS-S 
33ks ACIS-I 
merged ACIS-I 



1.0: 
1.0 



1.0+^ 



Q+fl:i 
"-0.1 
1 o+<'-2 



1.0 



5.4;f;l 
4.0:i;3 

□+0.9 
-0.7 

5 2+'-'' 
-'■^-1.4 
4,1+1.1 



3.9+ 



3.7 



-0.9 
+8.0 



-6.6 
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4.5+5-^ 
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54 
78 



83.0 
109.5 
130.9 

83.7 
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132.5 



Chandra ULX spectra 



BPL, two Hnes 


ACIS-S 
33ks ACIS-I 


0.3+** 

6+^'-^ 

05^:^ 

-Q-2 


-4.5+0 
-4.2+} 




merged ACIS-I 


-3.0+0 


BPL, no lines 


ACIS-S 


3+<'-' 


-4.5! 




33ks ACIS-I 




-2.6! 




merged ACIS-I 


5+<'-2 


-3.1!| 



7+0.8 

1.1+^-^ 
1.0!^'^ 

Q7+8:^ 

1 9+8i 

1+8:^ 

-0.3 



I 2+6-6 

1 7+^-^ 
'■'-1.7 

0+^-2 
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37 




4+^-6 

o.o!6^ 



5 
23 




48.9 
71.5 

103.3 
49.0 
71.5 
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XMM-Newton nuclear region spectra including point sources and the ULX 

BPL, three lines XMM-short 0.7!}; ] , 2.0!0 [ 12.6!<J f 6.6+11, 23 15.0!|^ 58 6.9+11 28 1024.9 1555 

XMM-long 1.6!};-],2.6!5]-^ 24.7!|'^ 1.3!|-5 9 8.1!]-^ 70 0.2!l-^ 2 1194.6 1573 



Note. — All errors are 90% confidence for 1 interesting parameter. Column 1: Abbreviations for the continuum model used, and whether the 6.4 and 
6.7 keV lines were fit for. Thermal Bremsstrahlung: BR; Power Law: PL; Broken power law: BPL; Bremsstrahlung with additional power law component 
representing unresolved X-ray binary emission: BR+PL; Power law with additional power law component representing unresolved X-ray binary emission; 
PL-l-PL; The break energy for the BPL model is not tabulated, but ranged from 5.6 to 6.3 keV in the best fit models to both the Chandra an d XM M-Newton 
spectra. The power law component used in the BR+PL and PL+PL models to represent the unresolved X-ray binary emission is described in 8 14.21 Column 3: 
Temperature kT (keV) for bremsstrahlung models, or photon index for power law models. Two photon indexes are given for the broken power law models fit to 
the ULX. Column 4: Continuum model normalization factor. These values are numerically lO'* X K, where K is the XSPEC model normaUzation factors. For 
the bremsstrahlung models K = 3.02 X 10"'5/47rD2 x J n^nidV where ne and «; are the electron and ion number densities respectively (cm"'), D is the distance 
to the source in cm and V the volume of the emitting region. For the power law and broken power law models K is the photon flux per unit energy (units of 
photons/keV/cm^/s), evaluated at £ = 1 keV Columns 5, 7 and 9: Line normahzations, effectively total line flux in units of 10"* photons s"' cm"^. Columns 
6, 8 and 10; Line equivalent widths corresponding to the best-fit line and continuum normalizations, in units of eV. Column 11: Best-fit value. Column 12: 
Total number of degrees of freedom. 
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Fig. 5. — Background-subtracted ACIS spectra of the nuclear diffuse emission, tlie summed nuclear point sources (excluding the ULX), and the ULX itself as 
seen in the ACIS-S observation (a) and the merged ACIS-I observations (b). To improve the clarity of the figure the enor bars are not shown on the point source 
and ULX spectra. The primary differences between the ACIS-S and ACIS-I spectra are purely instrumental, and are due to the poorer low energy sensitivity and 
spectral resolution of the cosmic-ray-damaged ACIS-I detectors. 
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Fig. 6. — Nuclear diffuse (panels a - 1), point source (excluding the ULX, panels m - o) and ULX spectra (panels p - r) and best-fit models with residuals. The 
first column are ACIS-S spectra, the second column spectra from the 33ks ACIS-S observation and the third column spectra from the merged ACIS-I observation. 
Panels a — c and g — i (the first and third rows) are pure power law and bremsstrahlung models respectively, in which residuals in 6.2 - 6.8 keV energy range are 
clearly visible. Panels d - f and j - 1 (second and fourth rows) are power law with lines and bremsstrahlung with lines model fits to the diffuse spectra. The point 
soiffce spectra are shown with best-fit models of a power law plus lines in panels m - o, and the ULX spectra with a broken power law plus lines model in panels 
p - r. The emission lines in the point source and ULX models are not required statistically. 
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Fig. 7. — XMM-Newton spectra of the nuclear region of M82, along with best-fitting models. The shorter of the two observations is shown in panel a, and the 
longer in panel b. Data from the PN detector is shown as filled black circles, MOSl data as open triangles and MOS2 data as open squares. The solid, dashed 
and dotted lines are the best fit model components. 
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tinuum model fits to the diffuse emission spectra (i.e. with- 
out any line components) show strong residuals in the energy 
range 6.2 - 6.8 keV in the ACIS-S, 33ks AClS-1 and merged 
(48ks) AClS-1 spectra, and the addition of the two lines does 
a good job of removing these residuals. The addition of the 
Gaussian model components representing Fe Ka and Fe Hea 
emission improves the quality of fit in all the diffuse emission 
spectra, but not in the point source and ULX spectra. Quan- 
titatively the best-fit values are reduced by between 7 and 
18 for the diffuse spectra, but change by only A^- = -1.6 to 
H-0.2 in point source spectra, and do not change at all in the 
ULX spectra. 

A power law continuum model gives a marginally better 
fit to the diffuse emission continuum than a bremsstrahlung 
model in ACIS-S and AClS-1 spectra, but visual inspection 
of the best fit models shows that the differences between the 
two best-fit continuum models are very subtle. It is notewor- 
thy that the slopes of the all of the best fit power law mod- 
els to the diffuse emission (2.0 < F < 2.7) are significantly 
different from the slope of the resolved point source spectra 
(F = 1 .0 ± 0.2). This is further evidence that the bulk of the 
diffuse hard X-ray emission does not arise in an unresolved 
population of point sources that are similar to but fainter than 
the resolved point source population. The power law slopes 
fit for the ACIS-S and AClS-1 diffuse emission spectra differ 
at greater than 90% confidence. We will address this issue 
when we fit the more elaborate models to the diffuse emission 
in §1321 

CCD-resolution X-ray spectra that fit well with a contin- 
uum plus weak line feature are often also well-fit by a con- 
tinuum plus absorption edge model. This model is physically 
plausible in the case of AGN where there may be very large 
column densities of absorbing material along the line of sight 
to the point-like continuum source. In the case of M82 where 
the apparent line feature is seen in the diffuse X-ray emission 
that covers a region several hundred parsecs in size an implau- 
sible amount of iron would be required in an absorption edge 
model. We therefor do not consider this case further. 

4.1.1. The significance of the iron lines in the diffuse emission 

spectra 

These spectral fits demonstrate the presence of weak line 
features in the diffuse emission spectra, and very weak or 
negligible line emission in the point source and ULX spec- 
tra. That the ACIS-S and AClS-1 spectra independently yield 
similar best-fit spectral parameters suggests that these results 
are not a statistical fluke or calibration oddity. Nevertheless 
we would like to better quantify the true significance of the 
Fe line detections in the diffuse emission, and the lack of de- 
tection in the point source and ULX spectra. 

In X-ray astronomy it has been traditional to use an F-test 
to assess the significance o f the addition of a weak line com- 
ponent to a spectral fit, but lProtassov et all (l2002l) show that 
this is inappropriate. We therefor used Monte-Carlo models 
to assess the significance of the line features in the diffuse 
and ULX spectra, as well as the uncertainties in the diffuse 
emission line fluxes. 

We created 40000 realizations of pure-continuum spectra 
(no lines) using both the ACIS-S and merged ACIS-I dataset 
best-fit power-law models, and fit them in the same manner 
as the real data with power law plus line models. We then as- 
sessed how often the fit-derived line fluxes in the 6.4 keV, 6.7 
keV and summed 6.4 and 6.7 keV lines equaled or exceeded 
the observed values. We found that we can reject the null hy- 



pothesis (that the lines found in the real data are merely noise) 
at high confidence for all lines in both datasets. 

For the weaker of the line features seen in the diffuse emis- 
sion, the E=6.4 keV line, the null hypothesis can be rejected at 
98.6225% confidence for the ACIS-S observation alone, and 
at 99.9225% confidence for the merged AClS-1 observation 
considered on its own. The probability that a spurious E=6.4 
keV line feature would appear in both sets of observations if 
there were no line emission at £ = 6.4 and 6.69 keV is essen- 
tially negligible. 

The null hypothesis can be rejected with even more confi- 
dence for the 6.7 keV line feature, at 99.99985% confidence 
for the ACIS-S data considered alone. In none of the 40000 
simulation of the merged AClS-1 data did Poissonian noise 
create a spurious 6.7 keV line feature with a flux equal to or 
greater than that observed. 

The E = 6.7 keV Fe line normalization derived from the 
ACIS-S data is larger than that from the AClS-1 data by a 
factor of ~ 1 .6, marginally larger than the 68% confidence re- 
gions. The continuum normalizations also differ in the same 
sense but by a factor of 1.2. Accounting for the continuum 
contribution from unresolved point sources (discussed below) 
will reduce the differences in spectral shape and normaliza- 
tion between the ACIS-S and AClS-1 spectra. 

4.1.2. Is only one of the lines real? 

We also considered the possibility that only one of the lines 
is genuine, and that the other is noise. Considering the case 
where only the 6.4 keV line was assumed to be genuine, the 
number of times we found a spurious 6.7 keV line with the 
observed flux in 10000 simulated spectra was 1 for the ACIS- 
S data and for the longer merged AClS-1 observation, i.e. 
this hypothesis can be rejected at > 99.99% confidence for 
each of the datasets independently. 

The chance that we could mistakenly find a 6.4 keV line 
feature if there were in reality only a 6.7 keV line feature is 
not negligible for the ACIS-S data: this hypothesis can only be 
rejected at 51% confidence. In the higher S/N merged ACIS-I 
observation this hypothesis can be safely rejected with 99.8% 
confidence. 

Note that these probabilities assume that the lines are intrin- 
sically narrow and are at the expected energies of E=6.4 and 
6.69 keV. If we allow the central energy (and optionally the 
width) of the Gaussian component represent the Fe Hea line 
to vary the statistical significance of the 6.4 keV line feature is 
reduced. Table|5]shows the best-fit line energies and values 
for models fitting for the Fe Hea line energy compared to the 
default spectral fits described above. We also experimented 
with removing the 6.4 keV line so as to assess whether a sin- 
gle line at £ ^ 6.69 keV provides as good a fit as two lines at 
£ = 6.4 and £ = 6.69 keV 

For the ACIS diffuse emission spectra fitting for the Fe Hea 
line energy provides essentially negligible improvement in the 
fit statistic. Removing the 6.4 keV line leads to marginally 
worse fits (Ax^ ~ 2). However a single broader line at E ^ 
6.5 - 6.6 keV provides as good a fit as two narrow lines at 
E = 6.4 and 6.69 keV in both the ACIS-S and merged AClS-1 
diffuse emission spectra. The best-fitting line width of a ^ 
0.2 ± 0.2 keV is to poorly constrained to be used as a physical 
constraint on the validity of the model. 

Based on purely statistical arguments we can not distin- 
guish between a model for the diffuse nuclear emission with 
two narrow lines at £ = 6.4 and E = 6.69 keV and a model 
with a single broader line at £ ^ 6.6 keV. We conclude that 
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TABLE 4 

Hard X-ray fluxes and luminosities 



Dataset 


Continuum, 2 - 


- 8 keV) 


6.4 keV Fe Ka Une 


6.7 keV Fe Hea line 


6.9 keV Fe Lya Hne 




fx X lU 


Lx X lU 


fx X lO'-* 


Lx X 10-" 


fx X 10"* 


Lx X 10"" 


X 10'" 


Lx X 10"3^ 




ergs"' cm"^ 


erg s"' 


ergs"' cm"^ 


erg s"' 


ergs"' cm"^ 


erg s"' 


ergs"' cm"^ 


erg s"' 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


ACIS-S 
33ks ACIS-I 
Merged ACIS-I 
XMM-short 
XMM-long 
ULX high state" 


2.86+Q 

2 76+''-" 
z./o„ 12 

T 81+011 


4 28+'' '' 


5.67!|^^ 
5.82;3;3Q 
5.53 

6 77+2'26 
°- ' '-2.30 
1 a',+0.94 
-0.94 
5 45+2.45 


8.78+1^^ 
8.58 

10.50+35^ 


11 00+^ °2 

Q T)+3.36 

O. / J j-j 

i6.io:|« 
8.7o:|:i» 

7 07+2.78 
'■^'-2.76 


17.00+"8 
11 0+5'^ 
13.50+^g 
24.90:3|5 

13.50:|;^^ 
11.40!4-3l 


<9.24 
<7.15 
<6.58 

7 74+2.59 
'•'^-2.60 
< 3.00 
7 47+2.79 
'•^'-2.62 


< 14.3 

< 11.1 

< 10.2 

12 0+" "2 

<4.66 
,1 f;n+"33 



Note. — Fluxes and luminosities quoted refer to the diffuse emission alone for the Chandra ACIS data. The contribution of the unresolved 
point sources has been removed from these values. Fluxes derived from the XMM-Newton data are the total hne fluxes from the nuclear region, 
which are a combination of the point source and diffuse X-ray emission. No broad-band flux is quoted for the XMM-Newton data as we do 
not attempt to separate the diffuse the diffuse and point source contributions. The diffuse emission fluxes shown above have been corrected 
from those derived from the diffuse spectra to account for the diffuse flux in the point source and ULX spectra regions (an increase by a factor 
1.32 for the ACIS-S data and 1.20 for both ACIS-I datasets). All errors are 68.3% confidence in one interesting parameter. Upper limits are 
reported at 99% confidence. Note that the E=2 - 8 keV energy band continuum luminosities do not include the hne emission. 
" This is the excess line flux in the 2001 May 06 observation (XMM-short) above that seen in the 2004 April 21 observation (XMM-long). 
This excess is presumed to come from the ULX in the high flux state. 



TABLE 5 
Best fitting Fe Hea line energies 



Dataset 


6.4 keV line present? 


Fe Hea line 


Fit statistic 






Centroid (keV) Width (keV) 


Ax" Au 


(1) 


(2) 


(3) (4) 


(5) (6) 



ACIS-S yes 6.66!|']° 

no 6.665™-'* 

no 6.6l!||-]^ 

Merged ACIS-I yes 6.69 ±0.17 

no 6.64:'<' * 

no 6.54^™-]^ 

XMM-long yes 6.65 ±0.02 

no 6.65 ± 0.02 

no 6.65 ±0.02 

XMM-short yes 6.66 ± 0.03 

no 6.66:'<' ['3 

no 6.63+^^"! 



(f) 


-0.3 


-1 


(t) 


1.7 





15+0-19 
'^'■^^-0.15 


0.2 


-1 


(f) 


-0.1 


-1 


(f) 


2.6 





20+"-" 


0.1 


-1 


(t) 


-6.7 


-1 


(f) 


-5.4 





0.07 ± 0.03 


-9.2 


-1 


(t) 


-2.3 


-1 


(t) 


5.6 





23+0-"** 
■ -0.07 


-2.9 


-1 



Note. — The change in the fit statistic with respect to the best-fit model from Tab[3]if 
the energy and width of the E ~ 6.7 keV is fit for, and if the E = 6.4 keV hne is present 
or not. For the Chandra diffuse emission spectra the best-fit model used was the "PL+PL, 
three lines" model. For the XMM-Newton spectra of the entire nuclear region the best-fit 
model was the "BPL, three lines" model. Confidence regions correspond to 90% confidence 
for one interesting parameter Column 2: Records whether or not a narrow hne at £ = 6.4 
keV is included in the spectral model. Column 3: The energy of the Fe Hea Gaussian line 
centroid (keV). Column 4: The Gaussian width of the Fe Hea Gaussian line (keV). Note 
that this is cr, not the FWHM. Models where this width was fixed at the default value 0.001 
keV are denoted with "(f)"- Column 5: The change in the best-fit x^ value from the default 
model in Table |3] A negative Ax^ indicates a better fit. Column 6: The change in the 
number of degrees of freedom v from the default model. 



the evidence for the £ = 6.4 keV iron line in the ACIS diffuse 
emission spectra is of marginal statistical significance. 

4.1.3. Iron lines and the point sources 

Including line components in the fits to the summed point 
source and ULX spectra leads to negligible changes in the 
best-fit values. From the existing fits it is safe to assume 
that there is no significant iron line emission in the spectra 
of the point sources (excluding the ULX). In the case of the 
ULX the apparent continuum level in the E = 6 -1 keV en- 
ergy band is significantly higher than in the diffuse or point 
source spectra. Would such a high continuum make it dif- 
ficult to see weak iron line emission in the Chandra ACIS 



spectra, or even miss line fluxes as high as those found in the 
diffuse emission? For the ULX to be a significant contributor 
to the total hard X-ray iron line emission from M82 it would 
have to produce a line luminosity comparable to that from the 
diffuse emission. We created 10000 simulated spectra of the 
ULX, again for both ACIS-S and ACIS-I observations, using 
the best-fit broken power law continuum model along with 
iron lines equal in flux to those in the best-fit diffuse emission 
spectra, and then evaluated how often we would have fit line 
fluxes for the ULX as low we observed. For the ACIS-S this 
would only tend to happen 15.3, 2.5, and 1.0% of the time, 
for the 6.4 keV, 6.7 keV and summed line fluxes respectively. 
The equivalent numbers for the ACIS-I dataset are 2.4, 0.7 
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and 0.1%. 

We conclude that despite the relatively high continuum 
level in the Chandra ACIS spectra of the ULX we could de- 
tect iron lines of comparable flux to that in the diffuse emis- 
sion spectra with relative ease. That we did not in either set of 
Chandra observations indicates that the ULX was not a source 
of significant 6.4 keV or 6.7 keV line emission in September 
1999 or June 2002. 

4.2. More elaborate fits the Chandra diffuse emission spectra 

It is worthwhile to consider two further elaborations of the 
previously performed spectral fits: accounting for the contri- 
bution of unresolved point sources to the continuum spectral 
shape, and placing limits on the line flux of Hydrogen-like Fe 
at £ = 6.97 keV. 

Based on the diffuse fractions described in § 13.21 we know 
that a small but non-negligible fraction of the Chandra diffuse 
spectra is contributed by unresolved point sources (URPS). 
The fraction of the apparently diffuse count rate in the E=3.0 
- 9.9 keV energy band arising in these unresolved sources is 
^diff-urps/Z^diff = l-(flff/.fdm) - 0.23 for the ACIS-S obser- 
vation and 0.13 for the merged AClS-1 observation. The 
longer exposure time of the merged AClS-1 observation leads 
to a significant improvement in point source removal, and re- 
duction in spectral contamination by unresolved sources. The 
power law slope of the summed spectrum of the resolved point 
sources (F ^ 1.0) is significantly flatter that the fitted dif- 
fuse emission continua (F - 2.1 (ACIS-S) or 2.7 (ACIS-1)). 
If the unresolved sources have a similar spectrum to the re- 
solved point sources then accounting for the URPS contribu- 
tion should steepen the fitted continuum spectral shape of the 
true diffuse emission. 

We modeled the contribution of the URPS to the diffuse 
spectra as a power law component with the same spectral 
slope as the resolved point sources, and with a normalization 
that produces the expected count rate /?diff-urps = (/diff-/diff) x 
/?tot in the 3.0 - 9.9 keV energy band. Given the uncertainties 
in fj^ff and R^ot this normalization is accurate to ^5%. The 
power law slope and normalization of the URPS component 
are fixed and are not allowed to vary during the spectral fitting. 
A Gaussian line component at a fixed energy of E = 6.969 keV 
was added to represent the Fe Lya line. The results of these 
fits are given in Table |3] 

Both bremsstrahlung and power law models were used for 
the continuum shape of the true diffuse X-ray component, al- 
though we found that a power law continuum was still favored 
statistically over a bremsstrahlung model. The best fit nor- 
malization of the 6.9 keV line is often zero and in all cases 
statistically consistent with zero (although note that the two 
sided confidence regions reported by XSPEC and given in 
Table [3] do not strictly apply when the region includes zero. 
Rigorously derived upper limits on the 6.9 keV line flux are 
described below). 

Best-fit values are lower (i.e. better) than the simpler 
models previously discussed, but the improvement is not sta- 
tistically significant. Nevertheless this exercise is worthwhile 
in that including the expected contribution from the URPS 
leads to a closer match between the fitted ACIS-S and AClS-1 
diffuse component continuum shape (7 or kT). 

We used Monte-Carlo simulations (with 5000 realizations 
each for the ACIS-S, 33ks AClS-1 and merged AClS-1 spec- 
tra) to derive upper limits on the 6.9 keV line flux and 6.9/6.7 
line flux ratio. Simulations were used to convert the limits on 
the line flux ratio into limits on the plasma temperature, as- 



suming collisional ionization equilibrium (CIE). These limits 
are shown in Table|6] and the corresponding line luminosities 
are given in Table |4] 

4.3. XMM spectra 

Nuclear region spectra were extracted from both of the 
XMM-Newton observations within a radius of 28. "6 of a = 
09h55"51?9,(5 = +69°40'47."1 (J2000.0), the same region as 
used for the Chandra spectra. It is impossible to separate 
the nuclear diffuse emission from the point sources and ULX 
emission with XMM-Newton, so the spectra shown in Fig.|2] 
represent all the emission from this region. We will not at- 
tempt to spectrally-separate the different components of the 
diffuse hard X-ray continuum, but are interested primarily in 
the Fe line fluxes and energies: are they consistent with the 
Chandra ACIS spectra; and can the superior sensitivity of the 
XMM-Newton spectra be used to detect or place stronger lim- 
its on the E = 6.4 and E = 6.9 keV line emission? 

For each observation the EPIC PN, MOSl and MOS2 spec- 
tra were simultaneously fit with the same spectral model, al- 
lowing for minor differences between the PN and MOS cal- 
ibrations by fitting for a multiplicative term that is applied 
to the model normalizations (typically the MOS 1 and MOS2 
model normalizations must be increased by ~ 5 - 8% [MOS 1 ] 
and 3% [MOS2]). The spectral data are fit only in the 3.4 
-9.0 keV energy range, and the spectra are rebinned to attain 
> 10 counts per bin. 

We find that a broken power law provides the best represen- 
tation of the continuum spectral shape in both XMM-Newton 
datasets. The best fit spectral parameters are presented in Ta- 
ble [3] the 6.9 keV line flux and 6.9/6.7 keV line ratios in Ta- 
ble |6l and the line luminosities in Table |4] 

In the longer XMM-Newton observation (2004 April 21) the 
E = 6.7 and 6.9 keV iron line fluxes are very similar to those 
seen in the Chandra ACIS spectra, although the 6.4 keV line 
flux is significantly lower in the XMM-Newton spectrum. 

Fitting for the energy of the Fe Hea line in this dataset 
results in a significant improvement in the fit at a line cen- 
troid energy of E = 6.65 ±0.02 keV (see Table |5]). Remov- 
ing the weak E=6.4 keV Fe Ka line results in slightly worse 
fit, but as with the ACIS diffuse emission spectra we find 
that the evidence for the presence of £ = 6.4 keV emission 
is of marginal quality. Once again a broad line (Gaussian 
fj = 0.07 ± 0.03 keV, 90% confidence in one parameter) pro- 
vides a better fit than a narrow Fe Hea line. This, if not due 
to systematic effects in the detector, corresponds to a velocity 
FWHM - 7400 ±3200 km S-' or the thermal Une width for 
r ^ 7 X 10'" K iron ions. These values are large compared 
to the maximum bulk velocities and temperatures expected 
in thermaHzed SN ejecta of v - 3000 km s"' and T - 10^ K 
(Chevalier & Clegg 1985), leading us to doubt whether the 
broadening of the E ^6.65 keV line is real. Unfortunately 
the lack of a strong (and intrinsically narrow) 6.4 keV line 
prevents us from ass essing the level of systematic effects (for 
example as done in T anaka et alJl2000() . 

In the XMM-Newton observation of 2001 May 06 the 
iron line fluxes are considerably larger than seen in any of 
the ACIS spectra or the longer XMM-Newton observation. 
Furthermore there is a significant detection of the 6.9 keV 
line, and the overall hard X-ray continuum is several times 
brighter than at the epoch of the XMM-long observation. It 
is most likely that the ULX is responsible for this additional 
flux, given previous and more extreme hard X-ray variability 
dPtak & Griffiflisll 19991; iMatsumoto & Tsurull 19991) . 
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TABLE 6 

The 6.9 keV line flux and 6.9/6.7 keV line ratio. 



Dataset 


6.9 keV Hne norm 


6.9/6.7 keV line ratio 


kTciE 


(1) 


(2) 


(3) 


(4) 


ACIS-S 


< 6.27 


< 0.80 


< 1.1 X 10** 


33ks ACIS-I 


< 5.39 


<0.97 


< 1.2 X lO** 


Merged ACIS-I 


<4.91 


<0.72 


< 1.0 X lO** 


XMM-long 


< 2.69 


<0.33 


< 7.6 X 10' 


XMM-short 


6.9 ±0.23 


0.46 ±0.27 


(8.7+2-0) X 10' 


ULX high state" 


6.7 ±0.39 


1.03 ±0.66 


(1.3™'^) X 10^ 



Note. — Upper limits are reported at 99% confidence. Where two-sided con- 
fidence regions are reported they are 68.3% confidence for one interesting param- 
eter. Column 2: Upper limit or best-fit value for the 6.969 keV line normalization, 
in units of 10"'' photons s"' cm"^. Column 3: Upper limit or best-fit value for the 
ratio of the photon flux in the 6.9 keV (Fe Lya) to 6.7 keV (Fe Hea) lines. Col- 
umn 4: Upper limit or value for the plasma temperature in Kelvin, derived from 
the 6.9/6.7 keV line ratio, assuming collisional ionization equilibrium (CIE) and 
using the APEC plasma code. 

" This is the excess line flux in the 2001 May 06 observation (XMM-short) above 
that seen in the 2004 April 21 observation (XMM-long). This excess is presumed 
to come from the ULX in the high flux state. 



Assuming that the iron Hne flux seen in the longer XMM- 
Newton observation comes from a non-varying diffuse X-ray 
component we can calculate the line flux due to the ULX 
alone in the high state it was in on 2001 May 06. For com- 
pleteness we present these estimated ULX line fluxes and a 
temperature estimate from the 6.9/6.7 keV line ratio in Ta- 
bles |6] & m even though they do not pertain to the diffuse 
emission. 

Fitting for the Fe Hea line energy in this observation only 
leads to a marginal improvement in the fit (Table|5]l. The pres- 
ence of the 6.4 keV line is strongly favored, if Fe Hea is nar- 
row, but as with the other nuclear spectra we have discussed a 
single broad line at £ ^ 6.6 keV provides as good a fit as two 
narrow lines at £ = 6.4 keV and E = 6.69 keV. 

5. DISCUSSION 

The multiple Chandra and XMM-Newton observations of 
M82 provide a relatively self-consistent picture of the diffuse 
hard X-ray continuum and iron line emission. The emission is 
brightest within the spatial region associated with the ongoing 
starburst and supernova activity. Emission at £ 6.7 keV (the 
Fe XXV Hea line complex) is detected with high statistical 
significance in all the observations we have analyzed. We can 
only place upper limits on diffuse E ^ 6.9 keV emission (the 
Fe Lya). The statistical significance of diffuse E ~ 6.4 keV 
emission is marginal in both the diffuse spectra obtained with 
Chandra and the total nuclear spectra provide by the XMM- 
Newton observations. 

5.1. Iron line emission: previous observations and point 
sources 

These results demonstrate that the iron line emission in 
M82 comes from a diffuse source and not from compact ob- 
jects in three of the four epochs sampled by Chandra and 
XMM-Newton observations. In general the iron line emission 
from the summed point sources and/or the ULX appears to be 
of negligible intensity compared to the diffuse iron line emis- 
sion. 

Nevertheless the distinctly higher iron line intensity (in par- 
ticular at £■ = 6.9 keV) seen in the XMM-short observation 
compared to the relatively consistent iron line intensities seen 
the in the other Chandra and XMM-Newton observations in- 
dicate that compact objects can contribute at times to the iron 



line emission from M82. It is most likely that the additional 

iron line emission came from the ULX. 

The analysis presented by Strohmaver & Mushotzkvl 
( |2003|) . where they assumed that all the iron line emission 
came from the ULX, was based on this XMM-short observa- 
tion. With our results we can now separate remove the time- 
constant diffuse contribution to the total iron line emission 
in this observation (Table |4]l. The 6.4 and 6.7 keV iron line 
fluxes from ULX in this high iron-line-emitting state were 
comparable to the total diffuse iron line emission. At this 
epoch the 6.9 keV line emission from the ULX is of com- 
parable intensity to the 6.7 keV emission. 

K-shell iron line emission has previously been detected in 
Ginga, ASCA and BeppoSAX X-ray spectra of M82, although 
the uncertainties in the mean line energy and flux were large . 
Ginga observations of M82 (lOhashi etalJI 19901: lTsun]lll992h 
detected emission at £ ~ 6.7 keV with an equivalent width of 
170 ±60 eV. Ptaket al. (1997) reported a marginally signifi- 
cant (3<t) detection of a line at £ ^ 6.6 keV with an equivalent 
width of 1 00 eV based on ob servations with ASCA. Using the 
same data lTsuru et al.l (1 19971) concluded that evidence for the 
line was not significant, and gave an upper limit at 99% confi- 
dence on the equivalent width of any line in the energy range 
6.3 <E (keV) <7.1 of less than 100 eV. Later observations of 
M82 with ASCA found a variable X-ray source (most proba- 
bly the ULX) in an extremely luminous state, and in these ob- 
servations the spectra displayed weak iron line emission with 
an equivalen t width of 121+gj eV (Matsumoto & Tsuru 1999) 



ICappi et all (il999h fit for both Fe Ka and Fe Hea fines in the 
BeppoSAX observations of M82, of which only the latter was 
detected (at E = 6.6 ± 0.2 keV) with an equivalent width of 
60!^^ eV. 

These historical values are similar to the range we find 
in the two XMM-Newton observations, where the total iron 
Hne equivalent widths were 109 eV (XMM-short) and 81 eV 
(XMM-long), although this may be a coincidence. Indeed, 
based on our results we would expect the iron line equivalent 
width in M82 to vary, as we know that the total hard X-ray 
continuum is strongly influenced by the variability associated 
with point sources, and also that at some epochs (in one of the 
four observations we analyze) one or more point sources also 
contribute to the iron line emission. 
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ICappi et al] (1 1999b provide Fe Hea and Fe Ka line pho- 
ton fluxes and upper limits in their analysis of the Bep- 
poSAX observation of M82 which are consistent with the iron 
line normalizations we find (Table O. The most convincing 
case for E = 6.4 keV Fe Ka emission we have is in XMM- 
short observation. Here the observed line photon flux was 
(6.6 ± 3.7) X 10~* photons cm"^ s"' , consistent with the up- 
per limit of < 1.5 X 10"^ photons cm"^ s"' in the BeppoSAX 
observation. The Fe Hea line flux is not well constrained 
in the BeppoSAX data: (1.6;|:[ o) ^ ^^"^ photons cm"^ s"^ 
which is consistent with both the flux from the diffuse emis- 
sion alone (^ 6 - 8 x 10~* photons cm"^ s"') or the larger 
flux of (l.S^Q j) X 10"^ photons cm"^ s"' in the XMM-short 
observation. 

Even in Chandra and XMM-Newton data the iron lines do 
not immediately stand out, so it is unsurprising (if disappoint- 
ing) that the observations of M82 by previous generations of 
X-ray telescopes do not place tight constraints on the histori- 
cal iron line luminosity of this galaxy. 

5.2. Expectations for diffuse hard X-ray emission in M82 

The presence of diffuse hard X-ray emission in the center 
of M82 from a T ~ 10^ K plas ma has been expected sinc e 
the seminal theoretical paper of lChevalier&CleggI ([1981 . 
Various authors have also considered inverse Compton (IC) 
scattering of IR photons by relativistic electrons (responsible 
for the strong non-thermal radio emission from the center of 
M82 ) as a significant source of diffuse hard X- ray emission 
(e.g. ISchaaf et al.lfT989tlMoran & Lehnertill997h . 

5.2.1. 6.7 keV iron line emission 

Diffuse thermal X-ray emission is expecte d in the standard 
super wind model (I Chevaher & Cleg3 119851 : iHeckman et alj 
Il990l) . and as long as there is little mass-loading (the mix- 
ing of cool er ambient gas into the merged SN ejecta, see 
ISuchkov e t al. 1996) the temperature of the hot plasma in 
the starburst region will be high enough for hard X-ray 
bremsstrahlung emission and 6.7 keV iron line emission. 

Here we will address whether the observed 6.7 keV iron 
line luminosities imply plasma properties consistent with a 
standard Chevalier & Clegg superwind model. We defer a de- 
tailed investigation of what constraints these hard X-ray ob- 
servations of M82 can place on starburst superwind models to 
a companion paper. 

The highest 6.7 keV line luminosities we can associate 
with the diffuse emission is Lx.e.v = 1-70 x 10''**ergs"' (Ta- 
ble |4|i. The implied starburst region plasma emission integral 
is EI = fnennV = Lx,6.7/(ZFe,*A6.7[7;,ZFe,0]) where the hot 
plasma of temperature T^, iron abundance Zpe.* (relative to 
the Solar iron abundance) and emissivity at Solar abundance 
A6.7[r,ZFe.0] occupies a fraction / of the starburst region vol- 
ume y*. 

Supernova and stellar wind ejecta in a starburst are rela- 
tively enriched in iron with respect to the Solar value, and 
at the temperatures concerned it is likely that there is lit- 
tle or no depletion onto dust grains. To determine Zpe,* we 
used version 4 of the Starburst99 population synthesis code 
jLeitherer et al.lll99 9). We used a continuous star formation 
model, with a Salpeter IMF between stellar masses 1-100 
Mq, forming stars of Solar metallicity, and evaluated the out- 
put once the stellar populations are essentially in a steady state 
balancing death and replacement 30 Myr after the onset of 
star formation). For t hese parameters the i ron ab undance is 
Zpe,* = 5.2Zfe on the lAnders & Grevessdd 19891) abundance 



scale used to derive the X-ray emissivities. Choosing a differ- 
ent time, or using an instantaneous burst model will change 
this value by less than ±50%. 

We assume that all of this iron is in the gas phase. Any iron 
initially depleted onto dust grains will be thermally sputtered 
in a time scale net ^ 10"^ yr cm"-' for plasma temperatures T > 
10^ K, and will have reached an ioni zation state equivalent 
to collision al equilibrium in this time (itohlll989l : ISmith et al.l 
1996; Bork owski & Szvmkowiakll997h . We will demonstrate 
below than n^t > lO'^yrcm"-' in the hot plasma at the center 
of M82, justifying our assumption of negligible depletion of 
iron into dust grains. 

The 6.7 keV Fe line emissivity peaks at kT = 5.4 keV (T = 
6.3 X 10^ K) at a value of Ae.v = 9.75 x 10"^^ erg s"' cm^ per 
unit Solar iron abundance (from version 1.3.1 of the APEC 
plasma code, see Smith et al. 2001). We will assume that the 
plasma is at this temperature (which is consistent with the up- 
per limits given in Table |6]l, but will demonstrate below that 
adopting a different temperature leads to minor changes in the 
derived plasma properties. 

The region occupied by the radio SNRs in M82 is ^ 700 pc 
in diameter, with a projected thickness of ~ 100 pc, yielding 
a total starburst region volume of = 1 . 14 x lO^-'cm-'. 

Adopting the simplifications that «h ~ 0.83ne and that the 
total number density /itot = 1.92ne, then the E=6.7 keV line 
flux in the ACIS-S observation implies that 

n,«0.19.r''M 7 V) 

(1) 

where the starred subscript refers to the default values adopted 
above. The filling factor of the very hot plasma within the 
starburst region is expected to be of order unity (unlike the 
soft X-ray emitting plasma in the larger scale wind), so we 
will assume / = 1 for all further calculations. 

Using these adopted values the thermal pressure in the star- 
burst region is approximately ~ ntot?" ~ 2.3 x 10^ K 
cm""* . This is consistent with estimates of the thermal pres- 
sure in warm ionized gas in the starburst region of P/k^ 1 - 

3 X lO^Kcm-3 (IHeckman et al.l[T990l: iSmifli et alJl2006h . al- 
though approxim ately an order of magnitude la rger than the 
values derived by iForster S chreibe r et al.l (1200 ll) . Preserving 
the linearity of the FIR-radio correlation between normal and 
starburst galaxies implies that the magnetic field strength in 
the center of M82 is B ~ 200 - 400^G (See Figs. 1 and 

4 of [Thompson et all 120061 Thompson, private communica- 
tion). The resulting magnetic pressure Pg/^ = B^/^nk ^ 1 - 
4 X lO^Kcm"^ would be in equipartition with the thermal gas 
pressure. 

The total mass and thermal energy associated with the hot 
plasma within the starburst region is M^, ^ ntotpniufVi, ^ 
2.2 X IO^'^Mq and ^jh,* ^ 1 .5n,„,kT^fV-, - 5.5 x lO^'^erg 
(where p. « 0.61 is the mean mass per particle in the plasma). 

The predicted continuum X-ray luminosity in the E=2 - 8 
keV band associated with this plasma is 2.49 x lO^*** erg s"' , 
approximately a factor 20 less than the observed diffuse con- 
tinuum. We will address this issue in greater detail below. 

Can the starburst in M82 create and maintain such a 
plasma? The speed of sound in the plasma is Cj sa 
i-fkT/fimay/^ « 1.2 X 10^ km/s. To zeroth order this over- 
pressured gas will flow out of the starburst region in a 
time"^ Tioss ^ R*/cs, with a mass flow rate M ^ M^,c^/R^, ^ 

' In practice gas takes longer than the flow timescale ri^ss =R*/ Cs to flow 
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TABLE? 

Plasma properties inferred from diffuse emission luminosities 



Parameter Units Dataset 

ACIS-S ACIS-S Merged ACIS-I ACIS-S Merged ACIS-I 

Input luminosity 6.7keVLx 6.7keVLx 6.7k:eVLx Continuum Continuum 

(1) (2) (3) (4) (5) (6) (7) 



Assumed kT 


keV 


5.40 


3.25 


5.40 


3.80 


2.70 


Assumed T 


K 


6.3 X 10' 


3.8 X 10' 


6.3 X 10' 


4.4 X 10' 


3.1 X 10' 


Assumed Zpe 




5.2 


5.2 


5.2 


5.2 


5.2 


EI 


cm 


3.5 X 10*' 


4.5 X 10''' 


2.8 X 10*' 


8.1 X 10*2 


1.1 X 10*3 


He 




0.19 


0.22 


0.17 


0.93 


1.08 


P/k 


Kcm-3 


2.3 X 10' 


1.6 X 10' 


2.1 X 10' 


7.9 X 10' 


6.5 X 10' 


Cj 


kms"' 


1190 


922 


1190 


997 


840 


M 


Mq 


2.2 X 10' 


2.5 X 10' 


1.9 X 10^ 


1.0 X 10* 


1.2 X 10* 


Era 


erg 


5.5 X lO'^t 


3.7 X lO^"* 


4.9 X lO'"* 


1.9 X 10^5 


1.5 X 10^5 


£KE/vfooo 


erg 


2.2 X lO'-* 


2.4 X 10^" 


1.9 X lO'-* 


1.0 X 10^5 


1.2 X 10'^ 


^-2-8 


erg s"' 


P 2.49 X 1038 


P2.16X lO^s 


P 1.98 X 1038 


' 4.43 X 103' 


' 4.35 X 103' 


^-6.7 


erg s"' 


' 1.70X 1038 


' 1.70X lO^s 


' 1.35 X 1038 


P 3.47 X 103' 


P 3.28 X 103' 


^6.9 


erg s"' 


P 4.02 X lO^' 


P l.OOx 10^' 


P 3.20 X 103' 


P 3.13 X 1038 


P 1.05 X 1038 


Tcool 


yr 


3.8 X 10* 


3.0 X 10** 


4.2 X 108 


7.2 X 10' 


6.3 X 10' 


Tflow 


yr 


1.7 X 10^ 


2.2 X 10^ 


1.7 X 10^ 


2.1 X 10' 


2.4 X 10^ 




Mq yr-' 


5.7 X 10-* 


8.2 X lO-"* 


4.6 X 10-^ 


1.5 X 10-2 


1.9 X 10-2 


Mflow 


M© yr-' 


1.3 


1.1 


1.1 


5.0 


5.0 




erg S-' 


1.0 X 10*2 


5.3 X lO'" 


8.9 X 10"*' 


2.9 X lO''^ 


2.0 X 10^2 



Note. — Plasma properties for the hot gas in the central region of M82 derived from either iron line or continuum 
l umin osities, and with various assumed plasma temperatures. The derivation of these parameters is discussed in 
§ 15.21 The luminosities quoted are either input (i) from the values given in Table |4] or predicted (p) based on the 
plasma properties inferred based on the input luminosity and spectral parameters assumed. Unit volume filling factor 
is assumed, and the volume and effective radius of the starburst region are = 1.14 X 10*3 cm3 ^n^j _ 210 pc 
respectively. 



1.3Mq yr"' (where we have adopted 7?* = (3y*/47r)'/^ = 210 
pc). To maintain the hot plasma the starburst mass injec- 
tion rate should be greater than or equal to this value. The 
associated starburst mechanical energy injection rate must 
balance the advective energy loss rate of £ ~ Eju.^Cs/Ri, ~ 
l.Ox lO'^^ergs-i. 

The IR luminosity of M82 is Lir = 5.8 x 10'" L©, 
based on the latest re-calibrated 12 - lO Oum IRAS fluxes 
dSanders et al.ll2003i) . Using the formulae in fKennicutJ ( ll998h 
the star formation rate associated with this luminosity is 
MsFR = 3.9M0yr"' (for a Salpeter IMF between 1 and 
IOOM0). For this star formation rate the starburst99 model 
described above predicts energy and mass injection rates of 
EsB = 2.4 X lO'^^eergs-i and Msn = 0.93/3 M© yr"', where 
e is the thermalization efficiency and /3 the degree of mass- 
loading. 

These numbers are not very sensitive to the assumed tem- 
perature of the plasma. In Table [T] we show the predicted 
plasma properties assuming a temperature of kT = 3.25 keV, 
the mean of thermal bremsstrahlung fits to the hard X-ray 
continuum. Plasma properties derived from the 6.7 keV line 
luminosity in the merged AClS-1 observation are also very 
similar to the numbers described above. In all cases the pre- 
dicted bremsstrahlung continuum luminosity is significantly 
less than that of the observed diffuse hard X-ray continuum. 

These estimates demonstrate that the observed E = 6.7 keV 
diffuse iron line emission from M82 is consistent with the 
starburst superwind model as long at > 25% of the SN me- 
chanical energy is harnessed by the wind. A more detailed 
comparison between model predictions and these observa- 

out of the starburst region. In lChevalier & Clegel 's ll98^ analytical model gas 
velocities within the starburst region are a function of radius, but are always 
subsonic. 



tions will appear in a companion paper. 

5.2.2. The mystery of the diffuse continuum 

Within 500 pc of the center of M82 we estimate that the to- 
tal diffuse X-ray luminosity is in the range Lx w (4.4 ± 0.2) x 
10^'' ergs"' in the £ = 2 - 8 keV energy band (see Table |4|. 
This luminosity has been corrected to (a) account for the dif- 
fuse flux excluded from the spectrum as because it fell in the 
regions with detected point source emission, and (b) account 
for unresolved point source emission. 

This val ue is hig her than the value of Lx = 2.2 x 10^^ erg s"' 
quoted by iGrifiiths et akl (|2000), but this difference is to be 
expected. The effective area of Chandra is now known to be 
lower than the values given in the early calibrations used by 
[Griffiths et al.l (l2000l) . and consequently the true luminosity is 
higher. 

It is tempting to interpret the diffuse hard X-ray continuum 
as the thermal bremsstrahlung emission associated with the 
iron line emission. Power law models fit the spectra statis- 
tically better than a bremsstrahlung model (Table [3]), but in- 
spection of the fits and their residuals show that the practi- 
cal differences between non-thermal and thermal models are 
small (Fig.|6]l. 

This interpretation is problematic for several reasons. The 
most fundamental problem with this hypothesis is that the 
equivalent width of the iron line emission would then im- 
ply a sub-Solar iron abundance if the plasma is in colli- 
sional ionization equilibrium. A secondary issue is that 
the plasma properties derived assuming the continuum is 
bremsstr ahlung emission be gin to strain the boundaries of 
a classic IChevaher & Clegg starburst model (even when ac- 
counting for emission from unresolved X-ray binaries, i.e. us- 
ing the BR-i-PL fits shown in Table O. 

The non-thermal hard X-ray continuum emission due to in- 



22 



Strickland & Heckman 



verse Compton (IC) scattering of IR photons off cosmic ray 
electrons has long been expected to be significant in star- 
burst regions dSchaaf et al.lll989l:lMoran & Lehnertlll997l;IOhl 
1200 lb . In the case of M82 we will show that IC X-ray emis- 
sion can not be a significant contributor to the diffuse hard 
X-ray continuum. 

5.2.3. Cotlisional equilibrium plasma models for the continuum 

The equivalent width of the diffuse E = 6.7 keV emission 
we measure in the central 500 pc of M82 is 0.24 - 0.32 keV, 
when corrected for unresolved point source emission. For a 
hot plasma with Solar iron abundance in collisional ionization 
equilibrium the 6.7 keV equivalent width lies in the range 1.6 

- 0.8 keV for temperatures in the range 7.5 < logT < 7.8, 
which would imply a gas phase iron abundance of Zpe ~ 0.15 

- 0.4ZFe,o. This is inconsistent with the super-Solar abun- 
dance expected from merged SN ejecta, implying either that 
the diffuse hard continuum is not bremsstrahlung from the 
iron-line emitting plasma, or that the plasma has had negli- 
gible enrichment beyond that of the ambient ISM (again this 
is under the assumption of collisional ionization equilibrium). 
As discussed above we do not expect depletion of iron onto 
dust to be significant given the densities, temperatures and 
flow time scales in the hot plasma. 

The metallicity of the stellar population and the warm 
ISM in M82 is surprisingly high given its moderate galactic 
mass; Stellar and H II region O and Ne abundances in M82 
are Solar or higher (lAchtermann & Lacylll995l: lOriglia et al.l 
12001: ISmiflietal.ll2006h . The are very few measurements 
of iron abun d ances , the only one we are aware of is that of 
lOrigUa et al.' ('2004') who derive a stellar iron abundance of 
Zpe = 0.46;!;q ^^Zpe,© from NIR spectroscopy"^. This is consis- 
tent with the largest values derived from the iron line equiva- 
lent width, so we will entertain the possibility that the diffuse 
continuum emission arises in a hot plasma. 

As can be seen in Table |7] the best fit bremsstrahlung tem- 
peratures (ACIS-S: kT = 3.8 keV; merged ACIS-I: kT = 2.7 
keV, see Table |3]l and luminosities (Table |4]l would indi- 
cate plasma volume emission measures well over one order 
of magnitude greater than those derived from the iron lines 
alone. Other derived plasma properties (e.g. mass, energy) 
are larger than those derived from the iron lines by a smaller 
factor, but still the energy flow rates equal or exceed the lev- 
els that can be supplied by the starburst. The implied thermal 
pressures are now higher than independent estimates of the 
pressure in the starburst region. These discrepancies alone 
suggest that the diffuse hard continuum is not predominantly 
thermal bremsstrahlung. 

Returning to the issue of the metal abundance it would be 
difficult for SNe to heat an ambient plasma to > 3 keV 
without mixing in metal enriched ejecta. To obtain a final iron 
abundance of Zpe ~ O.4ZFe,0 (i.e. very close to the ambient 
value) and hide any signature of enrichment by Zpe ^ 5ZFe,o 
ejecta would require a large mass ratio of ambient gas to 
ejecta, Mamb/Mej > 10. But if such heavy mass loading oc- 
curred the maximum gas temperature would be < 1 keV, 

We note that Origli a et al.1 120041) derived an iron abundance of Zpc = 
0.43g Qg^Zpe Q from the short XMM-Newton observation, for which we have 
shown the ULX contributes significantly to the continuum and iron line emis- 
sion. lOriglia et al. attempted to correct for point source emission using the 
Chandra ACIS-I observations, but given the difference in ULX states be- 
tween the different observations this would not have worked. It is purely a 
coincidence that their iron abundance estimate is similar to that we estimate 
for the diffuse emission alone. 



inconsistent with the spectral fits to the continuum and also 
inconsistent with the detection of any E=6.7 keV line emis- 
sion. 

We therefor conclude that it is unlikely that both the iron 
line emission and the diffuse hard continuum arise in a hot 
plasma in collisional ionization equilibrium. 

5.2.4. Non-equilibrium plasma models for the continuum 

If the iron ion populations are significantly far from col- 
lisional ionization equilibrium values it might be possible to 
reconcile the iron line emission with the diffuse continuum. 
The best-fit Fe Hea line energies in Table |5] although for the 
most part statistically consistent with the expected value of 
E = 6.69 keV, lie in the range E ^ 6.54 to E = 6.66 keV which 
is suggestive of a non-equilibrium plasma. Comparing the ob- 
served iron line fluxes with those predicted from the contin- 
uum we see that a Fe XXV ion population fraction about 1.3 
dex below that expected for a CIE plasma would be required. 

Electron and ion kinetic temperatures are far from equi- 
librium i n supernova rem n ants with high velocity sh ocks 
(see e.g. ' Vink et al.1 llOOl iRavmond & Korreckl I2005D . In 
a collision-less shock each particle species s thermalizes to 
a temperature 7^ = 2.4 x 10 (m5/mp)(v/100kms"')^, so we 
would expect iron to initially have a temperature significantly 
higher than the electron or proton temperature, before relax- 
ing into equilibrium via Coulomb collisions or on a faster 
timescale by plasma instabilities. 

For a electron number density of = 0.93 cm"^ and tem- 
perature T = 4.4 X 10^ K (values appropriate if the continuum 
is bremsstrahlung, see Table |7]i, the prot on-proton Cou lomb 
relaxation timescale is Tpp ~ 3.2 x 10^ yr (ISpitzeiill962l) . The 
electron-electron relaxation timescale is shorter by a factor 
(nip/meyi^ ~ 43. Assuming that the iron temperature was 
initially nife/mp times higher than that of the protons then 
the Fe-proton relaxation timescale is still short: Tfe-p ^ 400 
years. These time scales are short compared to the 10^ - 10^ 
year time scales significant within the starburst region, sug- 
gesting that the kinetic temperatures of the ions and electrons 
should be very similar Tfe ~ Tp sa T^. 

A given ion will approach CIE on a timescale 

[Wf (C, + a,)]"^ (2) 

where C, is the collisional ionization rate coefficient from 
charge state / to state /+ 1, and a,- is the total recombination 
rate coefficient out of state / into state /- 1 (Liedahl 1999). 
We calculated these rate coefficients fo r Fe XXV using the 
atomic dat a and methods given in in ,Amaud & Rothenflugl 
(Il985h and lArnaud & Raymond! (1 19921) . although we used a 
more accurate approximation for the first exponential inte- 
gral (Equations 5.1.53 and 5.1.56 in Abramowitz & Stegu^ 
11964 than the one given in lArnaud & Rothenflugl (I1985i) . For 
r = 4.4 X 10^ K we find UeTeq = 1.7 X lO^yrcm'^ These 
rate coefficients vary slowly with temperature, so hot iron 
in M82's starburst region should be close to CIE as long at 
>0.1cm-^ 

All of these estimates suggest that non-ionization equilib- 
rium conditions and order of magnitude changes in the Fe 
XXV ion fraction are unlikely in the iron-line-emitting plasma 
in M82. We can not explain why the best-fit Fe Hea line ener- 
gies are systematically lower than the collisional equilibrium 
value of 6.69 keV unless it is a statistical fluke, as systemat- 
ics would have to effect both the Chandra and XMM-Newton 
data. We are left with the conclusion that the diffuse hard X- 
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ray continuum is not thermal bremsstrahlung associated with 
the iron line emission. 

5.2.5. The continuum as inverse Compton emission 

The simplest non-thermal process capable of generating 
diffuse hard X-ray emission in a starburst is inverse Comp- 
ton scattering of IR photons by cosmic ray electrons. 

However, we do not believe that the dominant process cre- 
ating the observed diffuse continuum i n the center of M82 
can be IC X-ray emission. iKlein et aTl (Il988l) show that the 
radio spectrum of M82 has a spectral index of a ^ -0.3 at 
frequencies below 1 GHz, and -0.7 above 1 GHz. Assuming 
that mean wavelength of the IR photons is 80/im then IC X- 
rays in the 2-8 keV energy band are created by scattering off 
cosmic ray electrons with 7 in the range 360 < 7 < 720. Us- 
ingthe 1.49 GHz flux density of 7.66 Jy and size of 36" x 18" 
dCondon et al.ll9 90). the resulting minimum energy magnetic 
field strength is Bmin ^ 55^G (assuming the ratio of energy 
in protons to electrons is 100). This implies that the electron 
population responsible for the synchrotron radiation at ~ 20 
- 80 MHz is the source of any 2-8 keV IC X-ray emission. 
The resulting IC X-ray spectrum should have a photon index 
r = l-a~ 1.3. The Chandra observations constrain the X- 
ray continuum photon index F = 2.5 - 3.0 (±0.5 at 90% confi- 
dence), which is inconsistent with the F = 1 .3 slope that would 
be the signature of IC X-ray emission. 

We can also estimate the IC X-ray luminosity in the hard 
X-ray band given that the ratio of the IC to radio synchrotron 
luminosities is is equal to the ratio of photon to magnetic field 
density, Lxjc/^sync = Uph/Ub- To first order the IR pho- 
ton energy density t/pH ^ L/R/(47rc/?J) = 1.4x 10"^ erg cm"-' . 
The (minimum) energy density in the magnetic field is con- 
siderably smaller, Ub = -S~/87r = 1 .2 x 10"'" erg cm"-' . There 
is some uncertainty in what value to adopt for the radio syn- 
chrotron luminosity, given that the flatter spectral slope at 
frequencies below 1 GHz suggests absorption. Neverthe- 
less the emission is not optically thick, as the observed slope 
of a K, -0.3 differs from the a = +0.5 of an optically thick 
source. With no correction for absorption t he rad io spec- 
tral energy distribution given by iKlein et aH (Il988l) implies 
LsYNC = 2.0 X 10" ergs"! in the 20 - 80 MHz band. If 
we extrapolate from fluxes above 1 GHz assuming the in- 
trinsic spectral slope is a = -0.7 then we obtain Lsync - 
8.8 X 10^^ ergs"'. The predicted IC X-ray luminosities in 
the E=2 - 8 keV energy band are Lie = 2.3 x lO^'' erg s"' and 
L/p = 1 .0 X 10^^ erg s"' . These values are 4 to 20 times smaller 
than the observed diffuse X-ray continuum, and cement our 
conclusion that IC X-ray emission is not responsible for the 
di ffuse hard X-ray continuum in M82. 

[Thompson et al.l (l2006l) argue that the linearity of the FIR- 
radio correlation in galaxies indicates that the minimum en- 
ergy magnetic field strength underestimates the true magnetic 
field strength in starbursts. If this is true then our estimates of 
the expected IC X-ray luminosity are too high. 

The physical origin of the observed X-ray continuum in 
M82 is not clear to us at the present time. 

5.3. Similarities to the center of the Milky Way? 

Another peculiarity of the diffuse hard X-ray emission in 
M82 is that emission from nearly neutral iron (the E ^ 6A 
keV line emission) is not predicted in the standard models for 
X-ray emission from a starburst region. If the emission at 6.4 
keV in M82 is real and not a statistical artifact, then the X-ray 



emission from center of M82 does bear some resemblance to 
the diffuse hard X-ray emission seen in the Milky Way Galac- 
tic Ridge and Galactic Center. 

The Galactic Ridge X-ray Emission (GRX E, see e.g. 
[Warwick et al.ll980t lYamauchi & Kovama 1993) IS an appar- 
ently diffuse source of both hard X-ray continuum radiation 
and Fe line emission with a vertical scale height of ~ 100 
pc is found in the inner |Z| < 40° of the Milky Way. This 
emission has not been resolved into individual point sources 
in deep Chandra observa tions and h ence might be truly dif- 
fuse in origin (see e.g. Ebisawa et a l. 2001, 2005), although 
iRevnivtsev et al.l (l2006l) present a strong argument in favor of 
unresolved point sources as the origin of the GRXE. 

The GRXE has a total £ = 2 - 10 k eV energy band lu- 
minos ity of Lx = ( 1-2) X 10^** ergs"' ( Kovama etal] 1 1 9861: 
Tanurnaet al.lll999|) and a total iron line luminosity of ~ 9 x 
10^'' erg s"' (lYamauchi & Kovamalll993b . - 30 % of which 
comes from within ^ 2° (^ 560 pc) of th e Galactic Center 
(iKoyama et al.lll989t lYamauchi et al.ll990l) . These values are 
'-^ 5% of both the diffuse hard X-ray continuum and iron line 
luminosities we derive for M82. 

Line emission from £ ~ 6.4 keV (Fe Ka), E ^ 6.1 keV 
(Fe Hea) and E ^ 6.9 keV (Fe Lya) line is are often 
seen, although the relative line intensities vary with posi- 
tion (Kovama etal."1996'; ' Yusef-Zadeh et al.ll2002l: iPark et al.1 
2004; Muno et al. 2004). In fliis respect the GRXE and the 
emission from the Galactic Center (GC) differ from the nu- 
clear X-ray emission from M82, where the E = 6.9 keV line 
is not clearly detected. 

The spectrum of the hard X-ray continuum near the GC can 
be fit by either kT keV thermal models or F ~ 1 .9 power 
law models (IMuno et al |2004|) . This is a harder spectrum than 
the M82 hard diffuse continuum (A;r ~ 3 - 4 keV, F = 2.5 - 
3.0). 

Models which attempt to explain the GR and GC contin- 
uum emission as thermal bremsstrahlung from a. kT ~ 10 keV 
hot plasma are not currently favored (iTanaka et al.li2000l) . A 
plasma this hot is not gravitationally bound by the Galaxy, and 
given the large Galactic volume the GRXE appears to occupy 
any thermal model must to have a large mechanical energy 
input to balance the advective losses. Dogiel et al. (2002) es- 
timate, by assuming the continuum is thermal bremsstrahlung 
emission, that the mechanical power required in a thermal 
bremsstrahlung model for the GRXE isE ^ 10"*^ erg s"' . This 
exceeds the mechanical energy input from all supernovae in 
the Galaxy even assuming high thermalization efficiency e. 

It is also known that the intensity of E ^ 6.A keV 
line emission is hig hest in the vic inity of molecular 
clouds ( Kovam a et alj Il9 96; Murakamfitan |2000l 120011; 
lYusef-Zadeh et al.ll2002l;"1Park et al.ii2004 ^. The X-ray emis- 
sion in these regions is well fit by a reflection and fluores- 
cent iron line emission model, although none of the observed 
X-ray sources is currently bright enoug h to be the original 
source of the reflecte d X-ray emission ( Kovama et al.|[l996t 
'Murakami et al.''2000). Diffuse hard X-ray emission from the 
Sagittarius B2 cloud is well explained by irradiation from the 
direction of the GC by a transient source of X-ray luminosity 
Lx.2-1 3 X 10^'^ erg s"' (£ = 2 - 10 keV) (iMurakami et all 

llooih . 

Other models for the diffuse iron line emission from 
the GC include excitation of neutral iron by imp act ion- 
ization by low energy C osmic Ray electrons (IValinia et al.l 
i2000; Yusef-Zadeh et al.ll2002) or charge exchange processes 



24 



Strickland & Heckman 



dWargelin et al.l2005h . Another alternative is non-equilibrium 
emission excited in a ^ .6 keV plasma by quasi-thermal 
electrons ^Masai et al."2002'). The non-therma l mechanisms 
of Masai et al. (2002) and War gelin et"!!] (l2005h also produce 
E = 6.7 and 6.9 keV iron emission, and so might explain all 
the observed iron line emission. These exotic models appear 
plausible on energetic grounds, but more a set of rigorous pre- 
dictions and observations will be required to test these hy- 
potheses. 

So there are puzzling similarities and differences between 
the diffuse iron line and continuum emission in M82 and the 
larger scale GRXE in the Milky Way. High spectral resolution 
observations with the proposed calorimeters on Constellation- 
X are the only robust way to exclude the possibility that 
the iron line emission from M82 comes from an exotic non- 
thermal mechanism and is not wind-related thermal emission. 
It is clear that analogies between the normal spiral Milky Way 
and the extreme starburst superwind galaxy M82 must cease 
to be appropriate at some point. The existence of a large scale 
superwind in M82 directly requires that there should be a very 
hot plasma within the starburst region, unless our understand- 
ing of starburst-driven winds is deeply flawed. 

Deeper X-ray observations with the current generation of 
X-ray telescopes would be beneficial, in particular for better 
constraining the existence and fluxes of the weaker 6.4 and 6.9 
keV lines. We speculate that the diffuse Fe Ka emission from 
M82, if indeed present, might arise in a scaled-up version of 
the X-ray reflecti on nebulae seen near Galactic Center molec - 
ular clouds (e.g. iKoyama et aDll996l:lMurakami et al.ll2000l) . 
The most-promising radiation source for such a model would 
be the M82 ULX, which exhibits flares in X-ray luminosity 
up to Lx.2-10 6 X 10' * ° ergs"' , assuming isotropi c emission 
dPtak & GriffiflislfTggllMatsumoto & Tsuru|[T999h . 

With regard to the extended 6.7 keV iron line emission from 
M82 (but not the GC) we conclude that the most plausible ex- 
planation for it remains that it is thermal emission associated 
with wind fluid of merged SN ejecta in the starburst region. 

6. CONCLUSIONS 

We obtained new Chandra ACIS-S observations of M82, 
along with reanalyzing two earlier ACIS-I observations and 
two XMM-Newton observations from the public archive. 

Our aim was to re-inv estigate the diffu se hard X-ray emis- 
sion discovered by Grif fiths et al.l (l2000l) using the best cur- 
rent calibration and CTI correction methods, exploiting both 
the high spatial resolution of Chandra and the high sensitivity 
of the XMM-Newton EPIC detectors. 

With Chandra we can spatially separate the diffuse hard X- 
ray emission from resolved point source emission, and by us- 
ing the point source luminosity function we also correct for 
the contribution of lower luminosity unresolved point sources 
to the apparently diffuse emission. The XMM-Newton ob- 
servations lack the spatial resolution necessary for removing 
the point sources within the nuclear regions of M82, but pro- 
vide high S/N spectra of the summed diffuse and point source 
emission. 

Approximately 20 X-ray point sources down to a limit- 
ing luminosity of Lx ^ 4 x lO-'*' erg s"' (0.3-8.0 keV energy 
band) are detected within the central region of M82. We esti- 
mate that the ACIS-S observation resolves ~ 90% of all point 
source X-ray emission from M82 (^ 93% for the combined 
ACIS-I exposures). We remove and interpolate over the re- 
solved point sources to obtain estimates of the flux from truly 
diffuse soft and hard X-ray emission in M82. Correcting for 



unresolved point sources we estimate that 20 - 30% of the 
hard X-ray emission (3 < £ (keV) < 10) is truly diffuse, with 
a £■ = 2 - 8 keV luminosity of Lx,2-8kev ^ 4.4 x 10^^ erg s"' . 
The diffuse continuum is slightly better fit with power law 
models (best fit photon index F = 2.5 - 3.0) than with thermal 
bremsstrahlung emission (best fit fcr ^ 3 - 4 keV). 

The diffuse hard X-ray emission is brightest within the star- 
burst region, with a major axis extent well matched to the ex- 
tent of the region in which the compact radio sources (young 
SNRs) have been observed. The diffuse hard X-ray emission 
is relatively smoothly distributed within this region, sugges- 
tive of a larger volume filling factor than that of the soft X- 
ray emitting plasma. Strong, spatially-varying absorption is 
visible in the soft X-ray band within the plane of the galaxy, 
and this is likely to complicate X-ray-derived abundance esti- 
mates. 

The diffuse emission spectra from both ACIS-I and ACIS- 
S observations obtained within r = 500 pc of the nucleus of 
M82 definitely requires the presence of a £ ^ 6.7 keV line 
component (emission from highly ionized iron), and possibly 
the presence of £ = 6.4 keV Fe Ka emission (from neutral 
or nearly neutral iron). We used Monte Carlo simulations to 
assess the significance of these detections. The total nuclear 
region iron line fluxes in the 2004 April 21 XMM-Newton 
observation are consistent with those of the Chandra-derived 
diffuse component. We place upper limits on the diffuse E ^ 
6.9 keV Fe Lya line, which constrain the plasma temperature 
r < 8 X 10^ K. 

No statistically significant Fe emission is found in the 
summed X-ray spectra of the point-like X-ray sources or the 
ULX in the two epochs of Chandra observation. The total nu- 
clear region iron line fluxes in the 2001 May 6 XMM-Newton 
observation (previously studied in IStrohmaver & Mushotzkvl 
l2003h are significantly higher than the other Chandra and 
XMM-Newton observations. We attribute the excess iron line 
emission in this observation, where E = 6.4 keV and E = 6.9 
keV emission lines are also seen, to the Ultra-Luminous X- 
ray source in its high state. The iron line luminosity of the 
ULX in this flare state was approximately equal to the dif- 
fuse iron line luminosity of the starburst region. However, at 
most epochs of observation (when the ULX was not in a high 
luminosity state) the iron K-shell luminosity of M82 is domi- 
na ted by the diffuse component, co ntradicting the conclusions 
of lStrohmaver & Mushotzkvl2003l 

The luminosity associated with the diffuse iron lines is 
ix,6.7kev = (LI - 1.7) X 10^** ergs"' and Lx,6.4kev (2 - 
9) xlO-'^ergs"' (note that the 6.4 keV line emission is of 
marginal statistical significance). We consider the similar- 
ities and differences between the diffuse iron line emission 
in M82 and the Milky Way Galactic Ridge X-ray Emission, 
and conclude that a thermal model is the most plausible ex- 
planation for the 6.7 keV emission in M82. If we assume 
that the hot plasma generating the 6.7 keV iron line is en- 
riched with iron to the degree expected for SN ejecta in a star- 
burst (Zpe ^ 5ZFe,0) then the luminosity of the 6.7 keV iron 
line implies a gas pressure of P/k ~ 2 x lO^Kcm"-' , and re- 
quires mass and energy injection rates of M ^ 1.3 Mq yr"' and 
E ^ 1.0 X lO'^^ergs"' to maintain. These values are consis- 
tent with independent estimates for the starburst in M82, sup- 
porting our hypothesis that the diffuse 6.7 keV iron line emis- 
sion arises in the long-expected but previously unobserved 
wind fluid that drives the larger scale superwind in this galaxy. 

The nature and origin of the diffuse hard continuum is 
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harder to ascertain. It is difficult to convincingly interpret the 
continuum as the thermal bremsstrahlung associated with the 
6.7 keV iron line emission. For a plasma in coUisional ion- 
ization equilibrium this would require a low gas-phase iron 
abundance Zpe < O.4ZFe,0, much lower than the super-Solar 
abundance expected for supernova ejecta. The pressure and 
energy flow rates derived from the continuum also strain or 
exceed the standard starburst-driven wind model. Mass load- 
ing could not produce such a low iron abundance without 
cooling the gas to such a degree that no 6.7 keV line emis- 
sion would be produced. We explore the possibility of non- 
equilibrium ionization of the iron, but find that a variety of 
methods show that collisional equilibrium is expected in the 
iron-line-emitting plasma. Non-equilibrium conditions can 
not reconcile the continuum and iron Une emission. The steep 
slope of the continuum (which is quite different from the con- 
tinuum slope of the point sources, F 1) is strong evidence 



Abramowitz, M. & Stegun, I. A. 1964, Handbook of Mathematical Functions 
with Formulas, Graphs, and Mathematical Tables (New York: Dover) 

Achtermann, J. M. & Lacy, J. H. 1995, ApJ, 439, 163 

Adelberger, K. L., Steidel, C. C, Shapley, A. E., & Pettini, M. 2003, ApJ, 
584, 45 

Anders, E. & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 

Amaud, M. & Raymond, J. 1992, ApJ, 398, 394 

Amaud, M. & Rothenflug, R. 1985, A&AS, 60, 425 

Axon, D. J. & Taylor, K. 1978, Nature, 274, 37 

Bland, J. & TuUy, B. 1988, Nature, 334, 43 

Borkowski, K. J. & Szymkowiak, A. E. 1997, ApJ, 477, L49 

Bregman, J. N. & Irwin, J. A. 2002, ApJ, 565, L13 

Cappi, M., Persic, M., Bassani, L., Franceschini, A., Hunt, L. K., Molendi, 
S., Palazzi, E., Palumbo, G. G. C, RephaeU, Y., & Salucci, R 1999, A&A, 

350, 777 

Chevalier, R. A. & Clegg, A. W. 1985, Nature, 317, 44 

Chevalier, R. A. & Fransson, C. 2001, ApJ, 558, L27 

Colbert, E. J. M., Heckman, T. M., Ptak, A. F, Strickland, D. K., & Weaver, 

K. A. 2004, ApJ, 602, 231 
Collura, A., Reale, P., Schulman, E., & Bregman, J. N. 1994, ApJ, 420, L63 
Condon, J. J., Helou, G., Sanders, D. B., & Soifer, B. T. 1990, ApJS, 73, 359 
Cottrell, G. A. 1977, MNRAS, 178, 577 
Dahlem, M. 1997, PASP, 109, 1298 

Dahlem, M., Weaver, K. A., & Heckman, T. M. 1998, ApJS, 118, 401 
Dogiel, V. A., Inoue, H., Masai, K., Schonfelder, V., & Strong, A. W. 2002, 
ApJ, 581, 1061 

Ebisawa, K., Maeda, Y, Kaneda, H., & Yamauchi, S. 2001, Science, 293, 
1633 

Ebisawa, K. et al. 2005, ApJ, 635, 214 
Fabbiano, G. & Trinchieri, G. 1984, ApJ, 286, 491 

Forbes, D. A., Polehampton, E., Stevens, I. R., Brodie, J. P., & Ward, M. J. 

2000, MNRAS, 312, 689 

Forster Schreiber, N. M., Genzel, R., Lutz, D., Kunze, D., & Sternberg, A. 

2001, ApJ, 552, 544 

Freedman, W. L. et al. 1994, ApJ, 427, 628 

Frye, B., Broadhurst, T., & Benltez, N. 2002, ApJ, 568, 558 

Gehrels, N. 1986, ApJ, 303, 336 

GoUa, G., Allen, M. L., & Kronberg, R R 1996, ApJ, 473, 244 

Gonzalez Delgado, R. M., Leitherer, C, Heckman, T, Lowenthal, J. D., 

Ferguson, H. C, & Robert, C. 1998, ApJ, 495, 698 
Griffiths, R. E. et al. 2000, Science, 290, 1325 

Heckman, T. M. 1998, in ASP Conf. Ser 148, Origins, ed. C. E. Woodward, 

J. M. Shull, & H. A. Thronson (San Francisco: ASP), 127 
Heckman, T. M., Armus, L., & Miley, G. K. 1987, AJ, 93, 276 
— . 1990, ApJS, 74, 833 

Heckman, T. M., Lehnert, M. D., Strickland, D. K., & Armus, L. 2000, ApJS, 
129, 493 

Hickox, R. C. & Markevitch, M. 2006, ApJ, 645, 95 
Itoh, H. 1989, PASJ, 41,853 

Kaaret, P., Prestwich, A. H., Zezas, A., Murray, S. S., Kim, D.-W., Kilgard, 

R. E., Schlegel, E. M., & Ward, M. J. 2001, MNRAS, 321, L29 
Kennicutt, R. C. 1998, ARA&A, 36, 189 

Kilgard, R. E., Krauss, M. 1., Kaaret, R, Prestwich, A. H., & Ward, M. J. 

2003, Astronomische Nachrichten, 324, 121 
Klem, U., Wielebinski, R., & Morsi, H. W. 1988, A&A, 190, 41 



against Inverse Compton X-ray emission being the dominant 
source of the diffuse continuum. 

The logical next steps are to see what constraints these ob- 
servations place on theoretical models for superwinds (in par- 
ticular the efficiency of supernova heating the degree of mass 
loading), and test these models against observed 6.7 keV iron 
line fluxes for other starburst galaxies. 
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